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ABSTRACT 

The standard treatment of cooling in Cold Dark Matter hales assumes that all of the 
gas within a "cooling radius" cools and contracts monolithically to fuel galaxy forma- 
tion. Here we take into account the expectation that the hot gas in galactic halos is 
thermally unstable and prone to fragmentation during cooling and show that the impli- 
cations are more far-reaching than previously expected: allowing multi-phase cooling 
fundamentally alters expectations about gas infall in galactic halos and naturally gives 
rise to a characteristic upper-limit on the masses of galaxies, as observed. Specifically, 
we argue that cooling should proceed via the formation of high-density, ^ lO^K clouds, 
pressure-confined within a hot gas background. The background medium that emerges 
has a low density, and can survive as a hydrostatically stable corona with a long cool- 
ing time. The fraction of halo baryons contained in the residual hot core component 
grows with halo mass because the cooling density increases with gas temperature, and 
this leads to an upper-mass limit in quiescent, non-merged galaxies of ~ 10^^-Mq. 

In this scenario, galaxy formation is fueled by the infall of pressure-supported 
clouds. For Milky- Way-size systems, clouds of mass ~ 5 x IO^A/q that formed or 
merged within the last several Gyrs should still exist as a residual population in the 
halo, with a total mass in clouds of '^ 2 x 10^*^ Mq. The baryonic mass of the Milky 
Way galaxy is explained naturally in this model, and is a factor of two smaller than 
would result in the standard treatment without feedback. We expect clouds in galactic 
halos to be ~ Ikpc in size and to extend ^ ISOkpc from galactic centers. The predicted 
properties of Milky Way clouds match well the observed radial velocity distribution, 
angular sizes, column densities, and velocity widths of High Velocity Clouds around 
our Galaxy. The clouds we predict are also of the type needed to explain high-ion 
absorption systems at 2 < 1 , and the predicted covering factor around external galaxies 
is consistent with observations. 

Key words: Galaxy:formation — galaxies:formation — cooling flows — intergalactic 
medium — quasars: absorption lines 



1 INTRODUCTION 

Cooling and galaxy formation within dar k matter halos 
was fi rst ^ discussed in a modern context bv lWhite fc ReesI 
il978l) . who argued that gas cooling was a main driver be- 
hind the characteristic mass of galaxies. After halo collapse, 



^ Their ideas were based on th ose of iBinnevI llQTTT : 
iRees fc OstrikejiT97 7') and Sj^jm^ and were applied to CDM 



specifically by iBlumenthal et al.l il98'd) 



gas is assumed to shock-heat to the halo temperature, and 
to cool over a characteristic timescale 



kbT 

n,A(r) ' 



(1) 



which depends on the particle number density of the ion- 
ized gas rii, the gas temperature T, the Boltzma n con- 
stant fcj,, and the cooling function A(r). White fc FrenkI 
(|l991j) extended this approach in order to make predic- 
tions as a function of time and position in a halo. The 
framework assumes that all of the gas within a central, 
high-density "cooling radius" cools and falls in to fuel 
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galaxy assembly, while gas beyond this radius remains hot. 
This cooling-radius method for tracking gas cooling is cer- 
tainly a useful approach, and it has become the basis 
for gas accretion estimates in s emi-analytic galaxy forma- 
tion m odels (recently , jSomcrvill^^^L[j20WUBensonet^| 
[ 200,^: iHatton et al [2003; H ernauist fc SoringeJ l2003t 
iNaeashima fc Yoshii |2004). Hydro-dynamical simulations 
seem to verify thi s pictur e, at least roughly (Katz 1992; 



Thoul fc Weinberg 1995: Springel et al. 2001 : Yoshida et al.. 



1 



20021: iHellvet alJl2003ri . 

Implicit in this standard model is that all of the gas 
within the cooling radius cools and contracts monolithi- 
cally. It has been known for some time, however, that the 
hot gas associated with galaxies (and clusters ) should be 
thermally unstable and prone to fragme ntation (lFieldlll965i 
iFall fc Ree^ll985l: iMurrav fc Linlll990l) . At a given radius 
within a halo, the cooling time can increase for some gas 
and decrease for other gas as density and temperature dif- 
ferences become enhanced by the cooling process. The result 
is a fragmented distribution of cooled material, in the form of 
warm (~ 10*K) clouds, pressure-supported within a hot gas 
background. In this paper we argue that the consequences 
of including this expected ingredient may be far-reaching. 
The residual hot gas core has a low density, and thus can 
exist as a pressure-supported corona for a long time without 
cooling. The fraction of baryonic mass contained in the hot 
core component grows with halo velocity (or temperature) 
and we show below that this gives rise to a characteristic 
cooled, central galaxy mass of ~ 10^^ Mq in high-mass ha- 
los. Interestingly, this is roughly what is needed to explain 
the bright-end cutoff in the galaxy luminosity function. Sim- 
ilarly, including this multi-phase treatment can help explain 
the masses of Milky- Way type galaxies without the need for 
excessive feedback. 

In our picture, the gas supply into galaxies is gov- 
erned by the infall of warm clouds. We suggest that the 
cloud population will have a velocity dispersion similar to 
that of the host halo, and that clouds will fall in to feed 
galaxy formation only when cloud-cloud collisions or ram 
pressure forces rob them of angular momentum and kinetic 
energy. In Milky- Way-size halos, the residual population of 
clouds is expected to be substantial. The clouds occupy typ- 
ical galacto-centric radii of ~ lOOkpc, and can explain the 
High Velocity Cloud (HVC) population around the Milky 
Way and high-ion absorption systems in external galaxies. 
A characteristic cloud mass of ~ 5 x 10^ Mq matches most 
of the observed properties of H VCs (see J7|l . The same char- 
acteristic cloud mass is consistent with our theoretical ex- 
pectations (JHJ, aides in the understanding of absorption 
systems (3S|l, and helps explain the total mass of the Milky 
Way without any significant feedback {^^. Based on this ev- 
idence, we argue that there is direct observational support 
for the idea that gas fragmentation be included in models of 
CDM-based galaxy formation. 

In what follows we will use the Milky Way galaxy halo 
as a fiducial case for comparison. The properties of our 
" Milky Way" dark matter halo are adopted from results 
of iKlvpin et al. ( 2002) . The authors use a wide variety of 
Galactic data and a baryonic- infall calculation to deter- 
mine a best-fit halo mass of Mv — W^^Mq and an ini- 
tial halo maximum circular velocity of 'Knax ~ 163km s~^. 
Their Milky Way mass is motivated by the mass models 



of lDehnen fc BinnevI Jl998h . who obtain Mg = (4 - 6) x 
IO^Mq. 

The results of iKlvpin et all ll2002ll and 
IPehnen fc BinnevI lll99a) also provide a useful illus- 
tration of the "over-cooling problem" faced by the standard 
treatment of cooling in galaxy halos. In the standard model, 
the baryons that end up in the Galaxy are simply those that 
exist within the cooling radius: Mc ~ fcfbMv Here fc — 0.7 
is the fraction of baryonic mass within the cooling radius 
for Milky- Way type halos (see O and h = Qb/i^m =i 0.17 
is the cosmic baryon fraction JSpergel et al.M2003fl . Based 
on the numbers quoted above, the ratio of the "expected" 
cooled mass to the actual mass of the Galaxy, /g = Mq/Mc, 
is significantly less than unity: 



/g ==; 0.43 



Mg 



5 X 1010 Mf 





Mv 


-1 


\fbfc] 

L0.12J 


1012 Mo 



(2) 



If the standard cooling arguments are correct, then less than 
half of the baryons that have cooled onto the Galaxy still 
exist there today. If feedback is to explain this, it requires 
that the Galaxy lost half of its mass via strong winds with- 
out destroying the thin disk. In our picture, this difficult 
series of events is avoided because a large fraction of the 
mass within the cooling radius never fell in, but remains in 
the halo in the form of a warm/hot medium (see ijSJl. This 
effect becomes more important in high-mass halos because 
the density of the hot gas core (which scales like the cooling 
density) increases with halo temperature. 

We have made an effort to frame our results as an ex- 
tension of the standard treatment of cooling, and we give 
comparisons to the single-phase approach whenever possi- 
ble. Of course, gas cooling and accretion in galactic halos is 
more complicated than the static halo model we use. Indeed 
gas falling into halos may or may not be shock-heated effi- 
ciently, and some gas was likely accreted as cold material, 
either stripped from infalling satellites, or simpl y as "cold 
flows" jBirnboim fc Dekei20ollKeres et alj20CI4ll . Our goal 
is merely to extend the standard semi-analytic treatment to 
include an allowance for multi-phase gas during cooling and 
to work out the implications of this approach. We expect 
that no matter how gas clouds are accreted, our qualitative 
expectations will hold. Detailed, high-resolution hydrody- 
namic simulations will be needed to test these expectations. 
Unfortunately, as we discuss in Hill the numerical challenges 
facing such an endeavor may be significant. 

Before providing an outline of the paper, we men- 
tion tha t in a series of papers, D. Lin and collab- 
orators iJMurrav fc LinI 1199(1 Il992t iLin fc Murravl Il992l: 
Murrav ct al.' 1993: B urkert fc LinI 1200(1 iLin fc Murravl 
2000: .Murrav fc Litis :2004ll have explored the fate of warm 



clouds in a hot gas medium. The analys i s that follows builds 
on their work. iMo fc Miralda-Escudd (Il996ll consider the 
presence of pressure-supported clouds in dark matter halos 
in order to explain QSO absorption line systems. We make 
a similar connection in f|H| 

The next section contains a review of the properties of 
dark matter halos. In ij^i ^^ describe the standard treat- 
ment of radiative cooling in halos. ^Js^tsnds this model to 
include the formation of warm clouds within a background 
hot gas medium. In [JS] we explore cloud masses, bracket 
the allowed range, and discuss mass scales of interest. JH] 
is devoted to modeling galaxy fueling via cloud infall. In 
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Symbol 


Equation where first used 


Description 


fb 


0) 


The fraction of mass in the universe in the form of baxyons. 


M^,R^,Vv 


000 


Halo virial properties: mass, radius, and velocity. 


J^s 1 V^rnax ■, f-^v 


IHJ 


Halo parameters: NFW scale radius, maximum circular velocity, and NFW concentration. 


A(T),A, 


ca CHI 


The cooling function and its approximate metallicity scaling 


PcTlcRc 


in} o iTC} 


The cooling density, corresponding electron number density, and cooling radius 


f^i,fJ-e,Zg 


Gl} 031 


The mean gas mass per particle, mean gas mass per electron, and the gas metallicity. 


VT,Vd,VP 


^^ tm STI\ 


Ratios of average hot gas core temperature, density, and pressure to the values at the cooling 
radius. 


Pw,T^ 


03 113 


Warm cloud density and temperature. 


mcl.r-cl.Vci 


|23E3 03 


Cloud mass, cloud radius, and the typical cloud velocity. 


Mc,Mh,M,i,Mg 


i^ i^ El E5l 


The total mass in various phases: cooled material, hot core, clouds, central galaxy. 


Tram, Tcci T"in 


i43l l44t l47l 


Cloud population time scales: the ram-pressure time, the cloud-cloud collision time, the cloud 
infall time. 



Table 1. Frequently used symbols 



SQ and 35] we compare our expected cloud populations to 
HVC data and CIV absorption system observations, respec- 
tively. i|^ presents an explanation for the exponential cutoff 
in the bright-end of the galaxy luminosity function. Future 
directions and implications are discussed in HlOl and we sum- 
marize in m\ In what follows, we adopt a flat ACDM cos- 
mology, with parameters set by the best-flt WMAP values: 
h = 0.72, Qb h^ = 0.02 4, and Qmh^ = 0.14 iSoergel et alJ 
I2OO3I . see also lPrimacki (i2002) ) The implied fraction of mass 
in baryons is /;, =0.17. 



2 DARK MATTER HALOS 

A dark matter halo of a given mass Mv is cha racterized by a 
virial radius R^ . The spherical top hat model (iGunn fc Gotd 
Il972r) provides a reasonable estimate of the value of 7?v, 
which is set by the radius at which the average mass enclosed 
equals a characteristic virial density pv = A^pu- Here pu is 
the matter density of the universe and Av is a cosmology- 
dependent variable that can vary as a function of redshift. 
For our adopted ACDM cosmology Av — 360/(1 + z) when 
^^ 1 and Ay ~ 178 when z^ 1 (a more precise fit is given by 



iBrvan fc Norman loos') . With this definition the halo virial 
mass and radius are related via ^ 



Rv 



3A'U 
47rpv 



1/3 



K 






(3) 



■^ The virial radius and velocity scale with mass in the following 
way 

iJv ~ 206h- 

1/3 



Vv ~ 144km s 






The singular isothermal sphere (SIS) is a simple ap- 
proximation that is often adopted for the density profile of 
a dark matter halo: 



p{R) 



PvRv 

37?2 



* m 



47rG i?2 ■ 



(4) 



The SIS profile has a rotation curve that is flat as a function 
of 7?: V{R) — Vv — Vmax, where Vmax is defined to be the 
maximum rotation velocity of the halo. The temperature of 
a singular isothermal sphere is related to its velocity by 



t^I^^^^ = io'k( 



K„ 



■ykb 



\ 163km s 



(5) 



where Cg = Knax/\/2 is the sound speed of the gas, : 
proton mass, the polytropic index is 7 = 1 for an isothermal 
gas, and fii = 0.62 is the mean mass per particle (electrons 
and nucleons) of the ionized gas in units of the proton mass 
assuming a 30% mass fraction in Helium. 

While the SIS is convenient for illustrative purposes, a 
better fit to the res ults of cosmological N-body simulations 
is the NFW profile (iNavarro et alJll995l : iKlvpin et alJl2001^ 



p{R) 



psRt 



RiR + Rs 



(6) 



where pa is a characteristic density. The scale radius, Rs is 
often expressed in terms of the concentration Cv = Rv/Rs- 
Given a halo mass, the value of Cv sets the value of ps, 
and the profile is determined. Simulations show that the 
median value of Cv for a halo of mass Mv at redshift z is 
well-approximated by C^(M, z) ~ 9.6(M^/M,)-°-^^{l+z)-^ 
iJBullocket al.ll200ll ). Here M* ~ lO^M© is the charac- 
teristic mass for collapse ai z = for our cosmology (see 
iLacev fc Cola . 1993) . 

The maximum circular velocity for an NFW pro- 
file occurs at a radius -Rmax — 2.15-Rs, where V^ax = 
GM(7?inax)/i?max. For our adopted cosmology, a good fit 
to the virial velocity in terms of l^nax is 
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K ==; 0.468C 



(7) 

which is good to 1% for Vv between SOkm s~^ and 
1200kni s~^. We assume that in the absence of cooUng, the 
relationship between velocity and temperature for NFW ha- 
los foUows that for isothermal halos (equation |2J. 

Finally, N-body simulations show that the mass accre- 
tion history of a dark matter halo as a function of redshift z 
follows a rema rkably well-defined fu nction of the final con- 
centration C° dWechsler et alJl2002^ : 

(8) 



(km/s) 



M^{z) = Mv(0)exp 



C2 



Galaxy-mass halos have C^ ~ 13, and typically have ac- 
creted half of their mass by Zf ~ 1.1, corresponding to a 
lookback "formation ti me" of ^f ~ 8Gy r . 

As discussed in IWechsler et alJ (12002^ dark halos 
tend to grow qualitatively from the inside out (see also 
i Helrni et al.l2003l:IZentner fc Bulloc3l2003t iTasitsiomi et"ai] 
l2003ll . The central density remains roughly constant at the 
value set during an early, rapid accretion phase of halo 
buildup. Because of this, for a given halo, Knax is relatively 
constant, as a function of lookback time. Indeed for halos 
with Vmax — 175 — 250km s~^ at z = simulations show 
that Vmax stays approximately constant in the main progen- 
itor back to z ~ z/ (R. Wechsler, private communication). 
Motivated by these results, in what follows we will assume 
that the halo Vmax (and thus its temperature) will remain 
constant back to the time of formation. With these prop- 
erties of dark matter halos defined we can move on to the 
treatment of gas in a halo. 



3 HOT GAS AND COLLISIONAL COOLING 

In the standard picture of CDM-based galaxy formation, 
gas collapses with the dark matter, and subsequently 
shock heats to the t emperature of the virialized halo (see 
IWhite fc Frenkll99J) . The result is an extended halo of hot 
gas, which begins to cool over a characteristic timescale 
(equation 0, and provides the gas reservoir for galaxy 
formation. Recently it has been shown in hydrodynami- 
cal simulations that the situation is not as simple as this 
jBirnboim fc Dekel20o3lKeres et al.l2004ll . In low-mass ha- 



los, Vm 



< 



100km s~^ , the cooling time of the gas is so short 
that a shock can not be maintained. Thus gas accretes in a 
"cold flow". For halos with Vnax^ 100km s~^ simulations 
find that the gas does indeed shock heat to the temperature 
of the halo. In our multi-phase model of the gas, we can in- 
terpret these cold flows as gas that simply enters the halo as 
warm clouds (without first shock-heating and subsequently 
re-forming clouds). 

In what follows we work within the framework of the 
standard picture and investigate how the treatment of frag- 
mentary cloud cooling will change the results. We do so 
mainly in the spirit of direct comparison. Although, as we 
show below, most of the interesting ramifications of includ- 
ing multi-phase cooling occur in high-mass halos, where 
shock-heating is expected to occur, and the standard picture 
is roughly valid. As mentioned in the introduction, different 
expectations for the nature of gas accretion onto halos and 
the efficiency of shock heating will affect our results some- 
what, but the qualitative nature of our conclusions will not 



w) 100 



< 




1631 



10 
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Figure 1. The cooling function for five different metallicities as 
a function of halo temperature (lower scale) and corresponding 
halo maximum velocity (upper scale). The five metallicities, from 



top to bottom, are: Zg 



1.0,0.3,0.1,0.03 and 0.0 relative to 



solar. Also shown are simple power-law fitting functions described 
in Appendix ^ for each metallicity (dashed lines) . The power 
laws change slope at two characteristic temperatures, Tm and Tj, 
corresponding to when metal line and Bremsstrahlung cooling 
start to dominate, respectively. The fractional difference between 
the full function and the fit are shown in the bottom panel. 



change. For example, gas stripped from infalling halos will 
likely fragment into clouds during this process, as seen di- 
rectly in the the Magellanic Stream fe.g. lWeiner fc Williama 
llOOfi') and lead to a configuration similar to what we discuss 
below (although via a different chain of events). Of course, 
more direct modeling will be needed to test these expecta- 
tions in detail. 



3.1 The Initial Hot Gas Profile 

After halo formation, we assume that the hot gas obtains an 
extended density profile. Motivated by the non-radiative hy- 
drodynamic simulations summarized in iFrenk et alJ il999f) . 
we assume that a dark halo with an NFW profile and con- 
centration Cv will initially have hot gas that traces the DM 
at large radius, but that develops a thermal core at ~ 3-Rs/4: 



pKR) = 



R!po 



{R+^Rs){R + R. 



po = 



Mt 



47ri?f(?(Cv) 



(9) 



This profile gives a good fit to Fig. 12 in lFrenk et alJ il999li . 
The normalization, po, is set so that the total initial gas mass 
within the halo is Ali, — fhM^. The function g{x) describes 
the radial gas mass profile 

g{R/Rs) 



Ml{R)^Mb- 



where 



9{C.) 



g{x) = 9 In (1 + -a;) - 8 ln(l + x) 



4x 



(10) 



(11) 



l + x 

Note that the lCole et al.l ll2000l) group adopt a similar start- 
ing point for their hot gas profile (a non-singular isother- 
mal configuration). However, it has been common in other 
semi-analytic models to assume that the hot gas profile mir- 
rors that of the dissipationless dark matter halo profile. 
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Pg{R) = fbPdm(R)- We find that significant differences be- 
tween this approach (assuming an NFW profile) and the 
thermal-core assumption arise only for halos of cluster mass 
and above (where the cores become large). 



3.2 The Cooling Density 

Once the initial hot gas profile is in place, the time it takes 
for gas to cool is only dependent on its density and the 
rate of cooling, parameterized by the cooling function A(T). 
Since cooling is triggered by coUisional excitation, higher 
density material generally cools first. Given a time since the 
halo formed, tf, the "cooling density" is the characteristic 
density above which gas can cool: 



Pc = 



SfigiripkbT 
2fj,itfA{T, Zg 



(12) 



Assuming a 30% mass fraction in Helium"^ , the mean mass 
per particle in ionized gas is \ii = 0.62 and the mean mass 
per electron is /Xe = 1.18. In what follows we will associate 
the time f/ with the halo formation time, which can be es- 
timated using equation IHl 

The cooling function, A(r), can be calculated as 
a function of gas te mperature, T, and metallicity Zg 
JSutherland fc Dopital[l993l . We plot A(r) as a function 
of temperature for several different gas metallicities as solid 
lines in Figure The top axis shows the halo velocity that 
corresponds to temperature value shown on the bottom axis. 
The dashed lines show a series of power-law fitting functions 
presented in Appendix El 

Usefully, as long as the gas is only mildly enriched, 
Zg^ 0.1, then the cooling function in galaxy-size halos (with 
Tra < T < Tb) is well-described by a simple power-law form 



A(r, Zg) ~ 2.6 X 10" 



A, 



r 



W<^K 



3 --1 

cm ergs 



(13) 



The parameter Az is a constant that varies with the metal- 
licity of the gas. The lower limit on the range of validity 
corresponds to the temperature where metal line cooling 
starts to dominate, Tm — 1.5 x 10^, and the upper limit 
is set by the temperature when Bremsstrahlung becomes 
the dominate cooling process Tb ~ lO'^K -1-1.5 x ^9^''^10''K. 
We will typically concern ourselves with mildly enriched gas 
with Zg = 0.1, in which case Az = 1.0, and the above ex- 
pression is valid for halos with maximum velocities in the 
range ~ 60 — 300km s"'^. We comment that our metallic- 
ity choi ce is consis tent with metallicity estimates for some 
HVCs JTripp et al . 2003; Sembach et al. 2004). For other 
gas metallicities the values of A^ and the ranges of validity 
can be found in Table IXTl 

It will be useful to express the cooling density in 
terms of the corresponding electron number density ric = 



Pc/ifJ-e 



Adopting equation 1131 1 for A(r), we can derive 



a typical value for Wc using equation 1121 
6.1 X 10"^cm"^ T^iAzt^ 



(14) 



■^ A gas with 30% Helium by mass has 3 Helium atoms for every 
28 Hydrogen atoms. If the gas is fully ionized, there are roughly 
34 electrons for every 31 nuclei, with a mean mass per particle of 
/li ~ 8/13 = 0.62 and a mean mass per electron of /le — 20/17 = 
1.18. 
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Figure 2. The cooling density as a function of halo temperature 
for several choices of gas metallicity. The cooling density increases 
dramatically with halo temperature. Thus massive halos can re- 
tain a larger fraction of their baryons in the form of hot gas then 
low mass halos. 



We have scaled our results by the characteristic temperature 
and formation time of a Knax = 163km s~^, Mv — 10^^ Mq, 
Milky- Way type halo: 



106 K' 



*» = 8C^- (15) 

Here, 8 Gyr is the time since a halo of this temperature has 
grown by roughly a factor of 2 according to equation |H| 

We stress that the temperature scaling in equation 1141 
and all of the analytic expressions that follow are only ac- 
curate for galaxy size halos. A more accurate calculation, 
is shown in Figure |21 for several assumed gas metallicities. 
We have used the mass-doubling time in equation |H| to set 
the halo formation time at each temperature. This figure, 
and all of the figures to follow, rely on the tabulated cooling 
functions shown in Figure 

3.3 The cooHng radius 



IWhite fc Frenld il991f) applied the concept of the cooling 
density as a function of radius in a dark matter halo, and 
introduced the concept of a cooling radius as a method to 
track the amount of cold gas available to form stars. This 
method has been adopted by most subsequent semi-analytic 
models of galaxy formation and seems to do an adequate 
job of reproducing the results of hydrodynamic simulations. 
The cooling radius, _Rc, is defined as the radius where the 
cooling density matches the initial hot gas density: 



Pg(^c) 



(16) 



unless /9c < Pg(-Rv), in which case we will define Re = Rv 
For example, if the initial hot gas density follows that of a 
singular isothermal sphere then the cooling radius is 



Rr 



fbV,l 



^ttGpc 



217kpc r"'''^(AziJ 



sl/2 



(17) 



for pc > Pg(-Rv). 

While simple, the isothermal assumption is not a very 
good approximation for what we expect the initial hot gas 
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profile to look like (Figure|lJ and this can lead to large errors 
in the estimated size of Re (see upper left panel of Figure 
13 . If instead we use our adopted initial gas profile (equation 
1^ then the value of _Rc can be determined by solving a 
simple cubic equation. We find that an approximate fit to 
the solution for galaxy-size halos is 



i?c=il57kpc T-^'^iA^t^ 



Nl/3 



(18) 



This expression is a good fit for the exact cooling radius 
solution for halos with V^ax =— 120 — 400km s~^, as shown 
by the dotted line in the upper-left panel of Figure |3 (see 
below). For Vmax^ 120km s~^, the virial radius will be the 
relevant outer radius (long dash line in Fig. |3J : 



i?v==;253kpc T° 



(19) 



The above expressions are valid for 2 = in our cosmology. 
Figure 13 shows the exact solution to equations 1121 and 
ll6l for various assumptions for the initial hot gas profile, halo 
formation time, redshift, and gas metallicity (clockwise from 
upper left) as a function of halo Vmax. The solid line in each 
panel shows the derived cooling radius for our fiducial set of 
assumptions, and the long-dashed line in each panel shows 
the halo virial radius. As seen in the bottom left panel, gas 
metallicity is important in setting the cooling radius because 
metal rich gas can cool more efficiently (and therefore at 
lower density and larger radius) than metal poor gas (see 
Fig|2|l- Note in the upper left panel that there is a large 
difference between assuming an SIS initial gas profile (short 
dash) and our NFW-inspired assumption (solid). Not only 
do NFW halos fall off more quickly in density at large radius, 
but they have smaller virial radii for a fixed Vmax- 
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Figure 3. Cooling radius. Re, as a function of halo Vmax. In all 
panels, the solid line is Re calculated with our fiducial choices of 
initial hot profile (equation |^, metallicity {Zg = 0.1) and for- 
mation time (equation |3 at 2 = as a function of halo Vmax 
(assuming an NFW halo). Upper Left; Dependence of Re on the 
initial hot gas profile. The short-dashed line assumes an SIS gas 
distribution and halo profile. The "INT" label refers to our initial 
hot gas profile, and the dotted line is our power-law approxima- 
tion to this result feauation I18i . Note that an SIS halo has a 
larger virial radius at fixed Vmax. Upper Right: Dependence of 
Re on formation time. Lower Right: Dependence of Re on red- 
shift, for a fixed formation time (changing only the Vmax to -Rv 
relation). Lower Left: Dependence of i?c on metallicity. 



of hot gas that extends to the center of the halo and that 
can provide pressure support for the hot gas outside of Re- 



4 A TWO-PHASE MODEL OF COOLING 

In the standard prescription, the evolution of Re with time is 
used to evaluate the amount of gas available to form stars. 
All of the gas within the Re sphere cools into the central 
galaxy (over the halo formation time). The gas outside of 
this sphere is assumed to stay there, at the virial tempera- 
ture of the halo, tracing the background dark halo profile. 
As time goes on, the cooling radius grows, and so does the 
supply of cold, star-forming gas. In essence. Re is used as 
a book-keeping tool, since clearly this shell-like structure of 
hot gas represents an unphysical, hydro-dynamically unsta- 
ble configuration (at least if Re < Rv)- The physical situa- 
tion this approximation most closely mirrors is one in which 
all of the gas within the cooling radius cools and contracts 
monolithically over the cooling time, with hot gas from the 
outer regions moving in as a result. The implicit assump- 
tion is that the thermal instability inherent in the gas is 
unimportant in governing the gas infall within the cooling 
radius. 

A different, perhaps more physically-motivated picture 
arises by considering the two-phase nature of the gas. As 
mentioned, gas within Re is subject to the thermal insta- 
bility, and will tend to cool via cloud fragmentation. That 
is, not all of the gas within Re will cool, but rather a two- 
phase (warm/hot) medium will develop. Warm (~ lO^K) 
clouds will form and grow until the background density of 
hot gas is reduced to roughly pe . Thus there is always a core 



4.1 Cloud Formation and the Thermal Instability 



iFieldl ' 1965j first studied the thermal stability of astrophys- 
ical gases, which was extended to non-equilibrium systems 
bv iBalbusI (1198(11 . In the equilibrium case with no heating 
the instability criteria is 

ainA 



d\nT 



< 1. 



(20) 



Using the fitting formula described in Appendix ^ (equa- 
tion IA2I I. we see that hot gas should be unstable in all 
galaxy-size systems. Specifically, gas will tend to fragment in 
halos with temperatures above the metal-line cooling tem- 
perature Tm — 1.5 X lO^K. The range of instability extends 
below Tm if the gas metallicity is less than solar. 

The thermal instability leads to the rapid growth of 
perturbations and to the formation of warm gas frag- 
ments within the hot gas background (see H5.ll and also 
Murrav & Lin 1990). Perturbations can be either in the gas 
temperature or density and may be seeded by the accretion 
of substructure into the galaxy's halo. Most of the perturb- 
ing halos will have virial temperatures much below 10*K and 
therefore will not gravitationally bind the clouds that form. 
Of course, some of the most massive subhalos may drive per- 
turbations to become gravitationally bound to them; how- 
ever, as we will be focusing on a scenario with an order of 
magnitude more clouds than massive dark matter substruc- 
tures, we will assume that this is not the case for most of the 
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Figure 4. The solid line shows the initial gas profile given by 
equation ini for a halo with T = lO'^K (Vmax = 163km s'^). For 
comparison, an NFW profile (dot-dashed line) is shown normal- 
ized to the same total mass. The virial radius i?v, the NFW scale 
radius Rs, and the cooling radius Re are also marked. In the stan- 
dard cooling argument, all gas within the cooling radius contracts 
to form a central galaxy. In our two-phase model, there is a core 
of hot gas (dashed line) that extends to the center of the halo 
providing pressure support for the gas. In our scenario, the mass 
that cools is the integrated mass difference between the solid line 
and the dashed line. 



clouds. For a treatment of warm clouds in dark matter sub- 
structurc see Gnat & Stcrnbcrs (2004]) and lSternberg et alJ 
(|2002). 

The overdense, low-temperature regions that fragment 
will cool via atomic line cooling until they reach a temper- 
ature of T„ = 10* K and form clouds embedded within the 
hot, high-pressure background medium. Further cooling of 
the warm gas into cold (~ 300K) material likely will be pre- 
vented because of the presence of the extragalactic ionizing 
background, at least for the typical cloud densities that we 
derive below. 



4.2 Residual Hot Profile 

After the clouds form out of the original hot gas halo, a resid- 
ual hot gas component will be left. We can work out a model 
for this distribution by assuming that the gas returns to hy- 
drostatic equilibrium within the gravitational potential of 
its dark matter halo. If the residual hot gas does not radiate 
significantly then it will adjust to the pressure change adia- 
batically. This is roughly what is seen in the cores of the clus- 
ters simulated (Frenk et ah 1999) using non-radiative codes 
(note especially the high-resolution result of Bryant's code). 
Thus we assume that the gas is adiabatic within Re, with 
P oc pI^^ (as adopted by iMo fc Miralda-Escudelll99et) . Of 
course, it would be useful to test this assumption with more 
detailed multi-phase cooling simulations in the future. 

If we normalize by demanding that the hot gas reaches 
the cooling density at the cooling radius, and assume an 
NFW gravitational potential (neglecting the contribution of 
the baryons) then we find that the temperature and density 



Figure 5. The fraction of baryons in the halo that cool as a func- 
tion of halo maximum circular velocity in our model (solid line) . 
This fraction quickly falls from ~ 90% at 100 km 8"^ to ~ 2% at 
500 km s~^. Also shown is the same quantity computed using the 
standard, single phase cooling scenario (dashed line). The differ- 
ence arises from the additional hot gas core that develops in the 
multi-phase treatment. 



profiles of the residual hot gas halo follow (see Appendix 15} 



Phix) = pc 

Tkix) = T 



1 + — ln(l -f x) - ^ ln(l + Cc) 

1 + — ln(l + x)-^ ln(l + Cc) 
X Gc 



3/2 



(21) 



where the radius R is expressed as a; = R/ Rs and Cc = 
Rc/Rs- We have assumed that the hot gas temperature at 
7?c is equal to the halo temperature T, defined with respect 
to Vmax in equation|S] If Re < Rv, we assume that the profile 
outside of Re is isothermal.'' This solution is plotted for a 
halo with Vmax = 163km s~' in Figure |1| The presence of 
a hot gas core implies that at least some fraction of the gas 
within the cooling radius remains hot. 

For simplicity in the calculations that follow, we work 
under the approximation that the temperature, density, and 
pressure of the hot gas can be treated as constants as a 
function of radius within Re- As can be seen in Figure^] this 
is certainly a reasonable approximation for the density. In 
Appendixiniwe show that the volume-averaged temperature, 
density, and pressure of the hot gas within the cooling radius 
for expression 12H are given by 



Th = vtT, ph = r}dPc, Ph = VpPc, 

with rjT — 1.0, rjd — 1.35, and r]p 
values in our treatment below. 



(22) 
2.7. We adopt these 



* One may worry that this solution gives values of the hot gas 
density that are slightly higher than the cooling density at small 
radius. However, this is only true if this gas was sitting at this 
density for a time tf with the same temperature. As discussed 
in Appendix^ the gas in the core has likely fallen into the halo 
center more recently than tj, and was heated adiabatically as it 
fell. The higher gas temperature and the shorter time available 
for cooling will act to increase the cooling density of the central 
gas, and allow it to exist as hot material at a higher density than 
the global cooling density of the halo. 
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4.3 Cooling Efficiency 

Interestingly, including the simple expectation of a hot gas 
core changes the cooled-gas fraction in galaxy-size halos ap- 
preciably compared to the standard estimate. In the stan- 
dard model, all of the gas inside of the cooling radius is 
assumed to cool. For the initial gas profile given in equation 
|U]this mass is 

giCc) 



Mr^ = mI{<r^) 



Mb 



gic.) ' 



(23) 



where Mb — ftMv In our picture, Mh^ is divided between 
cooled gas and the hot gas corona. The total mass in cooled 
material, Mc, is always less than AIr^ because of the pres- 
ence of a hot core of mass Aih '■ 



Mc = M_ 



Ho 



Mh, Mh 



4 - r3 

-Trph-rtc- 



(24) 



Here ph oc pc is the average density of the residual hot gas 
profile within Re. 

The difference is illustrated explicitly in Figure |S] 
Shown is the fraction of baryons that have cooled in the halo 
for both the single phase (dashed) and multi-phase cooling 
(solid) as a function of halo Knax- For galaxy-size halos, 
the multi-phase treatment reduces the cold gas fraction by 
~ 40% compared to the standard case. This difference will 
be amplified if some fraction of the the gas that cools per- 
sists in the halo as warm clouds (see §|S||. The effect of the 
hot core is more important for high-mass halos because the 
cooling density increases with temperature (see Fig. |2J . For 
Vmax — 500km s~^ systems, the total amount of cooled gas 
is reduced by a factor of ~ 5 compared to the standard 
treatment. As discussed in !JHI this may have important im- 
plications for understanding the bright cutoff in the galaxy 
luminosity function. 



4.4 Cloud Size and Density 

For cloud masses of interest, self-gravity will not be im- 
portant in setting cloud sizes (see t|5.6t . Instead, pressure- 
confinement will set a typical cloud pressure and density. 
The implied density of a cloud is 



Pw = Pc- 



ypT 

T ' 



(25) 



where Tw is the temperature of the warm cloud. We assume 
that the clouds are roughly constant density, so that a cloud 
of mass rrici will have a characteristic radius 



'■cl = 



Snic 



4ttpu 



1/3 



: O.Skpc 



"^6 ^6 



'(A.i, 



Nl/3 



(26) 



where me = rrici / 10 Mq, and we have used r]p = 2.7 and 
Tw = 10^ K. 

At this stage, the cloud mass is the primary unknown 
parameter. We have normalized our cloud size using a char- 
acteristic mass rrici = 10^ Mq, and we argue below that this 
may be an suitable mass for a variety of reasons. Modeling 
the underlying mechanisms that determine cloud masses is 
beyond the scope of the current work. We will instead at- 
tempt to constrain the allowed parameter space of clouds 
using both theoretical limits and, later, observational hints. 
Specifically, in the next section we consider various processes 
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Figure 6. Properties of warm clouds as a function of Vmax. The 
upper shaded region shows the total mass within iJc as a function 
of the halo's maximum circular velocity. Clearly the cloud mass 
must be less than this. The solid line shows the Jeans mass of the 
cloud. The triple-dot-dashed line shows the minimum mass cloud 
that can form in the presence of conduction with /s = 0.2 (see 
il5.2l . Masses above the dot-dashed will survive conductive evap- 
oration over a time tj (again taking /s = 0.2). For masses below 
the short-dashed line, ram pressure drag will cause clouds to move 
more slowly than the halo velocity, making them unlikely candi- 
dates for High Velocity Clouds. The shaded region in the lower 
left shows where clouds will be destroyed by Kelvin-Helmholtz 
instability by having a cooling time longer than the cloud sound 
crossing time. The error bar shows the cloud masses of interest 
for Milky- Way size halos (Vmax = 163km s"'^) which help explain 
HVCs, high-ion absorption systems and the Milky- Way galaxy 
mass within our picture (see iQ JU and 86. 21 respectively'! . 



in the halo that will act to destroy clouds, and use these to 
set limits on viable cloud masses. 



5 CLOUD MASSES 

The physical processes that can act to limit the masses of 
warm clouds include conduction, evaporation, the Kelvin- 
Helmholtz instability, and the speed at which clouds can 
cool. Except for the last case, these processes impose a lower- 
mass limit on clouds that can survive. If the initial fluctua- 
tion distribution is power-law, it is perhaps reasonable to as- 
sume th at clouds will ten d to inhabit the lowest-mass regime 
allowed JLin fc Murravl fl992i. Another interesting mass is 
the Jeans mass, which does not necessarily affect cloud sur- 
vival, but may determine the mass scale above which star 
formation becomes efficient. Similarly, the relative impor- 
tance of pressure drag on a cloud compared to the gravita- 
tional force will vary as a function of mass and this will set 
a lower limit on mass scales of interest for HVCs. 

The discussion that follows is somewhat lengthy, and 
we provide a summary now, in conjunction with Figure |S1 
aimed at the reader who wishes to move beyond this section 
to the results. Figure |S| shows a space of cloud masses md 
versus halo Vmax. In order to guide the eye, we have placed 
an error bar on the figure to illustrate the cloud masses 
of interest for Milky- Way size halos (Vmax ~ 163km s~^) 
which help explain HVCs, high-ion absorption systems and 
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the Milky- Way galaxy mass within our picture (see ^ ^ 
and H6.2l respectivelv). The upper shaded region is excluded 
on physical grounds, as masses above its lower edge exceed 
the total baryonic mass within the halo's cooling radius, 
Mfic- The shaded region in the lower left corner, below the 
dotted line, is excluded because clouds in this mass- velocity 
regime will be destroyed by Kelvin-Helmholtz instabilities 
f H5.31 . The triple-dot-dashed line that runs just below the 
error bar is the characteristic cloud mass that arises if con- 
duction sets the cloud fragmentation scale in the initial hot 
gas halo ( ^5.'2t . The dot-dashed line shows the minimum 
cloud mass that could have survived evaporation within the 
hot gas halo ( H5.4t . The solid line is the dividing line between 
Jeans stable (below) and unstable clouds ( H5.6II . Finally, the 
dashed line shows the mass below which ram pressure drag 
will cause clouds to move at speeds below ~ Vinax, and thus 
be unlikely candidates for HVCs ( >I5.51 . The main conclu- 
sion here is that the cloud masses of observational interest 
are viable based on these considerations. 

The analytic expressions that follow were calculated as- 
suming the cooling curve power-law in equation 1131 and 
therefore are valid only for galaxy-size halos (Vmax — 60 — 
300km s^^). The lines in Figure |S| were determined using 
the Zg =0.1 cooling curve shown in Figure 



5.1 Cloud Formation: The Ability to Fragment 

As a result of the cooling instability discussed in H4.ll the 
contrast between temperature or density fluctuations in the 
initial hot halo will begin to grow as cooling proceeds. As 
slightly cooler regions begin to cool, they get denser and in 
turn, cool even more quickly, and this can lead to cloud 
formation. Specifically, if the over-cool region compresses 
more quickly than the background medium can cool, a sep- 
arate warm cloud wil l form within the hot gas background. 
iBurkert fc LinI i2000t) studied this process in some detail, 
and showed that warm, dense fragments will emerge in the 
hot medium as long as the density growth becomes nonlin- 
ear before the cooling becomes isochoric. The condition for 
cloud formation is that the sound-crossing time, t\ ~ A^/ch, 
along a perturbation of wavelength Xi, should be less than 
the characteristic cooling time for the halo, which by our 
definition of the cooling density equals tf. Here we have in- 
troduced Ch = Vmax/\/2 as the sound speed of the hot gas. 
Let us write the eventual cloud mass in terms of the ini- 
tial fluctuation size as rrici = 47r(Ai/2)^pc/3. The condition 
r\ < tf sets an upper limit on the cloud masses that will 
form 



important co mpared to cooling (or heat i ng) is know n as the 
Field length jMcKee fc Begelman|[l99ol:lFielJl965^ 



mci< 8.4 X W^^MqTJ^^A-^I 



(27) 



We conclude that all mass scales of interest should be able to 
form clouds before isochoric cooling occurs. Indeed, this up- 
per limit generally exceeds the total baryonic mass available 
within halos. 



5.2 Cloud Formation: The Conduction Limit 

A more interesting limit arises from considering conduction. 
If conduction is important in the hot gas halo, this can 
dampen temperature fluctuations and inhibit the formation 
of clouds. The length scale below which conduction will be 



Xi 



Tk{T) 



-,1/2 



niA(r) 



(28) 



where n is the conductivity of the gas. One can characterize 
the co nductivity as a fraction /s < 1 of the classical iSpitzeiJ 
il962l) conductivity: 



(29) 



_ L84xl0i°7^ _i _i 1 

K = /sKsp = /s ; — 7 — erg cm s K , 

In Ac 

where In Ac is the Coulomb logarithm, and we adopt 
In Ac ~ 35 as an approp riate value for the temp erature and 
density range of interest JCowie fc McKeg|l977^ . For an un- 
magnetized plasma, fs is unity and conduction is efficient. 
The presence of magnetic fields can make fs quite small, 
with ~ 0.001 if the fields ar e uniform or moderately tangled 
jChandran fc Cowlevll99a) . However lNaravan fc MedvedevI 
i2001r) have shown that /a ~ 0.1 in a medium where mag- 
netic field s are chaotic over a w ide range of length scales. The 
results of IZakamska fc Narava n (2003) imply that fs — 0.2 
can help solve the coo ling flow problem in clusters (see also 
iKim fc Naravanll2003fl . We will adopt fs = 0.2 as our fidu- 
cial value here. With this choice we find that the Field length 
of a hot gas at the cooling density is 

Xp ~ llkpc T^^^A'J\f^,^^\ 



(30) 



Scales smaller than Xp will tend to have a uniform temper- 
ature, and this implies a characteristic lower-limit on the 
mass: mJT; = 47r(AF/2)^pc/3. Using typical numbers we find 

m^, ^ 1.2 X lO^'Mo T^'/'AVhlf!^l (31) 

We plot the Field mass as a function of V^max for fs = 0.2 
as the triple-dot-dashed line in Figure |S| Note that m^i as 
defined above scales as /s . If conduction operates at ~ 
20% the Spitzer value, the characteristic mass scale is quite 
similar to our mass scale of interest. 

We point out that if Xf^ Rc, we expect cloud fragmen- 
tation to be stabilized completely. When this occurs, con- 
ductive heating from outside of Re can play an important 
role in setting the temperature structure of halos. This oc- 
curs when T^ 3.2 x 10^ K (1/max^ 920km s~^), or in massive 
cluster-size systems. (Note that the scaling in equation 1311 
is only valid for 60km s~^^ Knax^ 300km s"'^ because we 
have assumed a power-law form for the cooling function in 
its derivation.) 



5.3 Cloud Survival: 
Instability 



The Kelvin-Helmholtz 



Once clouds form they are subject to shearing stresses across 
their boundary as they travel through the hot medium. The 
flow can be subject to perturbations, and this is character- 
ized as the Kelvin-Helmholtz instability (KHI). 

The dominate destructive process is the "cham- 
pagne effect", which results from the development of a 
low-pr essur e, fast flow around the head of the movin, 
cloud dPor oshkcvic hfc Zeldovichlll98lt iMurrav et ai]ll99: 
IVietri etaLil997i'l . By Bernoulli's theorem, the pressure ex- 
erted by the fast wind at the head of the cloud is low, so the 
cloud's inner pressure can cause its material to be pushed 
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out from the top. lVietri et alj ilQQTT showed that the cham- 
pagne effect is stabihzed if the cooling time of a cloud is 
shorter than the sound crossing time of the cloud. That is, 
if the pressure waves inside the cloud are damped by ra- 
diative cooling before they can cross the cloud, then the 
inner part of the cloud cannot respond to produce the over- 
spilling. Note that this result holds even for clouds that are 
in thermal equilibrium with a background field, as assumed 
here. ^^^_^^^^ ^^_^ 

Based on the work of IVietri et alJ (|l993), we would like 
to compare the cooling time of our clouds to the sound cross- 
ing time. This will determine if they are stable against the 
KHI. The sound crossing time of the cloud is t^ = Tci/cw, 
where c„ — 11.5km s~^ is the speed of sound in the warm 
medium. Using equation 1261 for the cloud size we obtain 

r^'~6.7xl0Vs ml/%-\A,t,y/'. (32) 

We compare this to the the cooling time of a cloud of tem- 
perature T„ and density p„: 

cl 3/iempfci,Tw 1 c ^, in6 ^^-3 A j- (•qq\ 

'^'^ " o ^ \fT^ \ -^-^^^^ y^^ ^B ^^*8> (33) 

where we have used A{T^ = lO^A') = 4.9 x lO'^^cm^ergs"^ 
(see Figure 0. This allows us to define a characteristic 
KHI mass, above which clouds will be stable. By setting 
Tr^i {rrici) = Tc we obtain 



mS"'~lO.5M0 T-\A,t,Y 



(34) 



and note that the mass above which clouds are stable is a 
very strong function of the hot gas temperature. This is seen 
clearly by the shaded region in the lower left of Figure |S| As 
the host halo's temperature goes down, clouds become less 
dense (see equation 1251 , their cooling times increase, and 
they are more susceptible to the KHI. We see that cooling 
alone stabilizes most cloud masses of interest except in low- 
temperature halos.^ 



5.4 Cloud Survival: Conductive Evaporation 

Clouds may also be evaporated by conduction from the sur- 
rounding hot gas. The characteristic evaporation time scale 
is given by 

25fc(,mci 



lbTviiimpK(l jrci 
~ 16Gyr m^/^^-3/2(A,tJ-l/^ 

JCowie fc McKeelll977il where we have taken fs = 0.2. If 
we set this equal to the halo formation time (i.e. if require 
that clouds forming at f / have not evaporated by today) this 
gives us a lower bound on the cloud mass of 

3/2 
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25fc6 
3.5 X 10'' Mq Tl'^KV^tl. 



(36) 



This is shown as the dot-dashed line in Fig. El Cloud masses 
above this line will not evaporate over a time tf. 

^ Note that even in low-temperature halos, magnetic effec ts may 
stabi l ize clouds against t he KHI iChandrasckhar. .19611 iMiural 
ll984HMalagoli et al.lll996l) . 



5.5 Cloud Motion: Ambient Drag 

As the cloud moves through the hot gas halo at speed Vd, 
it will experien ce a "ram pressure" drag force that opposes 
its motion (e.g. lLandau fc Lifshitd^l959^ 

-Fram = -CdPhVclTTrd, (37) 

where Cd is the drag coefficient. Clouds reach terminal ve- 
locity, vt, when the gravitation force on the cloud is balanced 
by ram pressure force: 



2 _ ^ GM(D)mci _ 2 

^t — ^"Fv5 2 - ^ ^ ""^ma 

O^TrHjPhCd 



iDCiPh 



(38) 



For simplicity, we have assumed that the host halo is an 
isothermal sphere, with D the distance from the cloud to the 
halo center. If we assume that a typical distance is D ~ i?c 
then we obtain: 
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m^/V,^/-(A.tJ-/-, 



(39) 



where we have used rj^ — 1.35 and Cd = 1.0. Thus if clouds 
are sufficiently massive they can travel at a typical speed 
Vci — Vniax and not experience significant deceleration over a 
dynamical time. The limiting mass, below which Vd < Vmax 
is 



ram 
THd 



5.1 X lO-^Mo T-^'"{A,t,) 



a/2 



(40) 



This mass scale is plotted as the dashed line in Figure |S] 
Clouds very much smaller than this are physically viable, 
but their slow speeds would make them unlikely candidates 
for "high-velocity" clouds. We return to the effect that ram 
pressure drag will have on cloud motion in H6.ll 



5.6 Cloud Self-Gravity: Jeans Mass 

The Jeans mass estimates when a fiuid is unstable to self 

[ ravity. For a cloud confined by a pressure Ph it is given by 
Spitzeil[l97i) 



Mj = 



9ct 



5G'3/2pV2 • 



(41) 



Warm clouds with temperature Tw — 10* K have a sound 
speed Cw ~ 11.5km s~^, and the minimum cloud mass that 
is unstable to self gravity is 



■m-ci 



' ~ 7.2x 10**Afo T}^/^'^-^ ^"-'^ 



'(A.i. 



(42) 



The Jeans mass for a warm cloud as a function of halo max- 
imum circular velocity is plotted as the solid line in Fig. 

n 

Our expectation is that clouds will be less massive than 
the Jeans mass in galaxy-size halos and therefore their self 
gravity can be ignored. However, it is possible that in high- 
mass halos that clouds will tend to be more massive. This 
might happen, e.g., if conduction sets the cloud mass (equa- 
tion l31l and the triple-dot-dashed line in Fig.|njl. In this case, 
the clouds may be above the Jeans mass, self-gravitating, 
and perhaps form stars in cluster-size halos. These objects 
would likely be spheroidal systems, with rather low mass-to- 
light ratios. Cluster "galaxies" of this type would not contain 
dark matter. 
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6 GALAXY FORMATION VIA CLOUD INFALL 

The warm clouds are accelerated towards the center of the 
halo, but we do not expect that they will settle there im- 
mediately. Rather, clouds should have some distribution of 
energy and angular momentum that will have to be lost be- 
fore the clouds merge with the central galaxy. For example, 
clouds may take on an angular momentum distribution simi- 
lar to that of bulk- ave raged regions seen in dark matter halo 
N-body simulations ll BuUock et al.l 120011) or o f streaming 
motions in hydrodynamic simulations (Ivan den _ Bosch et al. 
l2003l) . Clouds stripped fror n merging halos are ex pected to 
have a similar distribution JMaller fc Dekeil2002ri . 

For concreteness, we will assume that the clouds take on 
an energy distribution with a characteristic velocity equal to 
the halo velocity, Vd — Vmax- Cloud-cloud collisions or ram 
pressure drag eventually lead to cloud orbital decay. In our 
picture, it is this cloud infall that sets the gas supply that 
governs the formation of the central galaxy. In the following 
subsections, we discuss expected infall times, explore the 
implied central galaxy mass and residual cloud population, 
and discuss the disruption of clouds as they approach the 
galaxy. 



6.1 Infall Times 

As clouds orbit within the hot gas halo they will experience 
ram pressure drag (equation I37t . The continuous drag will 
sap energy from the clouds, and this can eventually lead to 
orbital decay. The timescale for this to occur is given by 

2mci 



TrCdr^iPhVci 



2.6Gyr ml^'^T- 



1/2 



(A,t, 



a/3 



(43) 



if Vci = Vmax and Cd = 1- Thus for 10 M© clouds, only 
those that cooled out of the hot gas more than ~ 3Gyr ago 
would have begun to sink to the center of the halo via ram 
pressure effects. If the drag coefficient, Cd is less than 1, 
then this is a lower limit on the drag decay time.® We show 
below that the mass in clouds that we expect to have fallen 
in or formed in the halo since that time may be a rather 
large fraction of the baryonic content of the Galaxy. 

Cloud collisions will also be important in triggering 
cloud infall. We will assume, for simplicity, that after a col- 
lision, most of the cloud energy goes into heating the cloud 
material, and that this is quickly radiated away. The remain- 
ing, likely merged, system will have low kinetic energy and 
will quickly fall in to contribute to the central galaxy. Thus 
the cloud infall time will scale like the cloud-cloud collision 
time. 

The mean free time between cloud collisions can be writ- 
ten as Tec — (<?f>ciVciO"ci)~^, where the cloud cross section is 
(Tci ~ vrr^i and 4>ci is the number density of clouds. Note that 
(j}ci depends on the total mass of warm clouds Ma, and that 
this will change as clouds collide and merge or if new clouds 
form. If we assign cj>ci — 3Afci/(mci47r_Rc) we can write the 
cloud-cloud collision time as 

® We note that the Reynolds number for clouds in such a halo, 
Re~ fciVci/i' (where i/ is the viscosity), is expected to be quite 
high, Re~ 10^, if conduction sets the viscosity. In this case, the 
drag coefficient, Cd, is likely to be less than unity, even for super- 
sonic flow. 



3MciVcir2j 



2.4Gyr ml/'T^^''{A,t,)'^'A4^,lo.m 



Here we have used M2.10 = Mci/{2 x IO^^Mq) as the charac- 
teristic mass in clouds that we expect to exist in our fiducial 
halo. We explain this expectation in more detail in H6.2I We 
stress, however, that Md should vary as a function of halo 
mass and cloud mass, so their are additional dependencies 
in equation I44l that are hidden in this variable. 

When the density of clouds is high, Tcc is small, and 
clouds will quickly collide and sink to the center. As the to- 
tal number of clouds drops, the cloud infall rate will begin 
to drop as well. It is useful to consider the simple scenario 
where we start with a number density of clouds <j)o in a 
fixed volume. In this case the number density of clouds as 
a function of time obeys d(f>/dt — 4>{t)T~^ = (fi'^{ro(j)o)~^ , 
where tq is the mean free time initially. The solution is 
0(t) — (j)o/{l + t/To), so there will always be a residual cloud 
population in any halo, if cloud collisions set the infall rate. 

We mention that the shortest timescale over which 
clouds can fall in to the central galaxy is the free-fall time, 
Tff. If we estimate the cloud free-fall time from the cooling 
radius as th — -Rc/Knax then we obtain 



rff~0.94Gyr r-'/**(A,tJ 



Nl/3 



(45) 



This expression is accurate for halos with maximum circular 
velocities of 120km s"-^< Vmax^ 400km s"-^ (where _Rc < -Rv 
sets the lower limit and the breakdown in the scaling A(T) ex 
T~^ sets the upper limit). For all cases that we consider, the 
free-fall timescale is shorter than both Tcc and Tram. 



6.2 Central Galaxy Mass 

The total mass within the halo cooling radius, Mu^ , is di- 
vided between gas in the hot halo core, Mh, and gas that 
has cooled since the halo formation time Mc (equation I24II . 
The cooled gas mass is itself shared between warm clouds, 
Mci, and the central galaxy Mg. The mass budget is then 
described by 



Mc = Aid + Mg 



(46) 



We assume that cooling proceeds by the formation of 
clouds and that the infall of clouds leads to galaxy growth. 
The evolution of the total mass in clouds as a function of 
time can be modeled as a competition between cold mass 
accunrulation (as the cooling radius grows) and the rate of 
cloud "destruction" via infall onto the galaxy: 

dMci _ dMc dMg 
^ dt 



dt 

dMg 

dt 



dt 



(47) 



Mci 

Tin ' 



We have associated the cloud infall rate with the rate of 
galaxy growth, and set this equal to Md/nn- Here Tin is a 
characteristic cloud infall time, chosen to be the minimum 
of Tcc and Tram. For the fiducial halo and cloud mass dis- 
cussed in this section, Tcc < Tram, and cloud-cloud collisions 
dominate the infall. 

It is straightforward to solve this simple set of equations 
1461 I47II in order to evaluate the mass in each component. 
Once we choose a halo Vinax, the evolution of the hot gas 
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Figure 7. Evolution of eacli baryonic component as a function of 
lookback time. The long daslied line shows the the total baryonic 
mass within the cooling radius, which grows steadily as a func- 
tion of time. The central galaxy mass (solid line) also grows with 
time, as clouds continue to fall in. The total mass in the hot gas 
core (short dashed line) remains roughly constant as the cooling 
radius grows and cooling density drops. The total mass in clouds 
(dot-dashed line) remains nearly constant, as cloud formation and 
infall are balanced. 



core mass, Mh(t), is governed entirely by the evolution of 
the cooling radius and corresponding evolution in the cool- 
ing density (equations 1141 UHl 1241 . The evolution of Mc with 
time can similarly be determined by the evolution of Rc{t) 
(eauation ll8ll and Mh(t) (as just described). The other com- 
ponents may be tracked via equation 1471 once one adopts a 
cloud mass, rrici, and evaluates rin(i) using equation 1441 or 
US 

Figure |7| shows the resulting buildup in each mass com- 
ponent as a function of lookback time for our fiducial "Milky 
Way" halo of Vmax = 163km s~^, and a cloud mass of 
rrici = 5 X 10^ Mq. The top long-dash line shows the to- 
tal baryonic mass within the cooling radius, Mr^, and the 
lower set of dot-dashed, short-dashed, and solid lines show 
the galactic mass, hot core mass, and total cloud mass (Mg, 
Mh, and Md) respectively. The most striking result is that 
the final galaxy mass, Mg ~ 6 x 10^^ Mq, is roughly half 
of the mass it would have been had we adopted the stan- 
dard treatment, and allowed all of the mass within _Rc to 
contribute, Mn^ ~ 12 x 10^° Mq. The mass in clouds peaks 
rather early at Md ~ 2 x 10^° Mq, and then remains rela- 
tively constant as the cloud infall rate is matched by the rate 
of accumulation of cooled gas. The mass within the hot core 
also remains relatively constant as a function of lookback 
time Mh ~ 4 x 10^°Mq. This is because Mh oc ricRc, and 
growth in Re at late times is canceled out by the decrease 
in ric (this can be seen via inspection of the time scalings in 
equations ^] and 1181 . 

The final residual cloud mass expected in this scenario 
is straightforward to understand without having to solve 
the differential equations. Given a characteristic cloud infall 
time Tin, we expect clouds to remain in the halo as long as 
Tin is longer than the time since the clouds were formed (or 
accreted) . 

Consider a case where rin ~ Tcc and cloud-cloud colli- 



Figure 8. Shown as the solid line is the amount of mass that has 
cooled as a function of lookback time: Mace (i) = ^c (0) — Mc (t) . 
The dashed lines show the total mass in clouds that would have a 
collision time, Tec, that is equal to the lookback time plotted on 
the X-axis. We show the result for three cloud masses, (1, 5, 10) X 
10^ Mq. The point where a dashed line intersects the solid line 
gives the approximate mass in clouds that can survive until t = 
(today) . For a typical cloud mass m^i = 5x10^ Mq , the expected 
residual mass in clouds is Md ~ 2 x 10^^ Mq (as indicated by the 
dotted line), which is in good agreement with what is found by 
integrating the set of equations discussed in the text and shown 
in Fig.|7| 



sions dominate infall. The solid line in Figure |H| shows the 
amount of mass that has cooled since a lookback time t. 
Specifically, M^cc{t) = Mc(0) - Mc{t). We have used the 
same Knax = 163km s~^ halo discussed in conjunction with 
Figure [7| The three dashed lines in Figure |H| correspond to 
cloud-cloud collision times computed using three different 
cloud masses md ~ (1, 5, and 10) x W^Mq. In each case 
we show the amount of mass needed in clouds to get a Tcc 
equal to the time on the x-axis. Note that cloud-cloud col- 
lision times are short if the total mass in clouds is large. 
The point where the dashed lines cross the solid line corre- 
sponds to the mass in clouds that could have survived until 
the present day without merging into the galaxy. Thus, for 
the central line {rud — 5 x 10^ Mq) we expect a final mass 
in clouds of Md — 2 x lO^'^M© to have survived until the 
present day. ^From Figure |7| we see that this is very close 
to the mass in clouds derived from integrating equations 1461 
and l47ll . This simple method allows a quick way to estimate 
the final central galaxy mass: Mg(0) ~ Afc(O) — A/acc(T). 

We have applied this simple treatment in Figure |5| in 
order to estimate the range of cloud masses that may help 
explain the baryonic mass of the Milky Way. The shaded 
band show s the estimated range of M ilky Way galaxy masses 
derived bv lPehnen fc BinnevI lll99a) . The solid line shows 
the final galaxy mass, Mg, and the dashed line shows the 
residual total mass in clouds, Md, as a function of the in- 
dividual cloud mass, nid- Note that the Milky Way mass is 
matched well for md — 4 x 10® — 10* Mq without includ- 
ing any blow-out feedback. Note however, that even if cloud 
masses are small (rrid^ 10® M©) and fall in quickly to as- 
semble the Galaxy, explaining the mass of the Milky Way 
is much easier in this picture because of the substantial hot 
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Figure 9. Shown is the mass in the hot gas core (dashed), mass 
in warm clouds, with and without ram pressure (dot-dashed and 
dotted) and mass in the galaxy (solid) as a function of the warm 
cloud mass. The shaded region shows the range of masses for the 
Milky Way (^Dehncn & Binnev 1998). Cloud masses in the range 
(4 — 100) X 10^ Mq allow agreement with estimated Milky Way 
mass without including any blow-out feedback. 



gas core. The kinks in the lines occur at the cloud mass be- 
low which ram pressure dominates cloud infall. Cloud-cloud 
collisions are more important for massive clouds. The dotted 
line shows how the mass in clouds would change had we set 
the cloud drag constant to zero, so that ram pressure was 
unimportant in cloud evolution. 



6.3 Tidal Disruption 

As clouds approach the central galaxy they will experience 
strong tidal forces and may be destroyed. If clouds are bro- 
ken up before impacting the galaxy, this will prevent them 
from causing too much heating or damaging the disk. Tidally 
deformed clouds will tend to be quite large, and perhaps 
can be identified with the large HI complexes t hat are well 
known to exist in proximity to the Milky Way ( Tripp et alJ 
l2003l: IWakker fc van Woerdenlll99ll: iBlitz et alj lToog 



Clouds will be disrupted when the tidal force from the 
host potential overcomes the pressure confinement of the 
clouds. At a distance D from the host center, assuming that 
the host potential is roughly isothermal, the tidal force felt 
across a cloud of radius rd ^ D is approximated as 

GMdm(-D)mcirci 



F- 



D3 



This can be compared to the typical pressure force 
Fp ~ Tir^iTippcT. 



(48) 



(49) 



The radius R^ where the two forces are equal defines the 
surface of a sphere of disruption for the clouds. For an SIS 
density distribution M{R) = G^^V^^^R, the tidal sphere 
for destroying clouds can be derived analytically: 

^^(A.tJ^''^ 



R. 



13kpc 






(50) 



In the galaxy mass regime (100km s~^ < Vmax < 350km s~^) 
the value of Rd for a SIS and a NFW halo are very sim- 
ilar. We see that the sphere of disruption is larger than 



the size of the Milky Way's disk. It is therefore unlikely 
that the disk will be heated significantly from impacting 
clouds. Interestingly, this distance would be consistent with 
distance limits for many o f the large HVC com plexes (e.g. 
Ivan Woerden et aljlT999bllat IWakker et alJl200j) . 



7 RESIDUAL CLOUDS AS HVCS 

Seen in 21 cm HI emission. High Velocity Cl ouds (HVCs) 
have been studied for more than four decades JMuller et alJ 
1963). Interpreting the observed properties of HVCs in terms 
of physical parameters requires knowing the radial distance, 
D, from the Sun. This is the major observational unknown 
that has fueled the debate over their origin since their dis- 
covery. 

Models for HVCs range fro m condensed "Galactic 
Fountain" gas at D ~ 5kpc dShaoiro fc Fieldl 119761 : 
iBregmara Il980l) . to large, extra-galac tic objects assoc i ated 
with the Local Group D ~ IMpc iVerschuuilll969l : lArd 
F985; Blitz et al.' 1999; "Blitj l2002l : ISternberg et al.l l2002l: 
iMalonev fc Putman 2003) . In our picture, the HVCs are 
"circumgalactic" , within the cooling radius of the Galaxy, 
and bear resemblance to the D ~ lOOkpc population sug- 
gested by Oort (196^. "^ 

We expect that most of the mass in each cloud is in 
the form of ionized hydrogen at a temperature T ~ 10^7^. 
Clouds of this kind would have a velocity distribution fuU- 
width-half-max (FWHM) of Aw ~ 27km s"^ The hne width 
distrib ution of Compact HVCs studied by Ide Heii et alJ 
J2002h has a median FWHM of Av = 25km s'^. The agree- 
ment with predicted and observed line widths is encourag- 
ing. 

The recent HIPASS survey cataloged HVCs over the en- 
tire southern sky IJPutman et al. 2002). They find that the 
radial velocity distribution of their clouds is narrow when 
plotted with respect to the Galactic Standard of Rest, with 
(Tr — 115km s~^ (it is ar — 185km s^^ for the Local Stan- 
dard of Rest). It peaks* near ~ 0km s~^. As discussed in 
the introduction, completely disjoint dynamical models of 
the Milky Way lead us to choose a fiducial "Milky Way" 
dark halo with Knax = 163km s~^, and therefore a veloc- 
ity dispersion very close to the GSR distribution of HVCs: 
ar — Knax/\/2 = 115km s~^. Since the residual clouds in 
our scenario should roughly take on the dark halo's velocity 
distribution, we expect them to match the observed HVC 
distribution quite well. 

The HIPASS HVC population has a characteristic peak 

~^ an d a characteris- 

"aPEooa) . 



itmanet alJ 1200211 . The 
expected hydrogen space density for an individual cloud in 
our model is 



HI column density of Nhi — 10 
tic angular size of S ~ 0.5deg'^ JPr 



^ We mention that a fragmentary, pressure-supported HVC pop- 
ulation similar to the one we suggest might arise with a Local 
Group barycenter, as long as Andromeda and the Galaxy share a 
common hot gas halo (L. Blitz, private communication). Here we 
will focus on the Galaxy as an isolated halo as we have throughout 
this work. 

** The HVC velocity distribution set in the Local Group Stan- 
dard of Rest peaks near ~ —75km s~^, and has about the same 
dispersion as the Galactic Standard of Rest distribution. 
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Uh = 



pwfh 



1.4 X 10" cm"^ T^(Azts 



(51) 



where we assume that the mass fraction in Hydrogen is 
fH = 0.7. The column density in HI through a cloud can 
be estimated via Nhi = ir^inH^Hi, with eni the fraction 
of neutral hydrogen. With eni =0.1 IjMalonev fc PutmarJ 
[2003) we obtain 



Nhi ^ 6.8 x 10 ci 
If we demand Nhi 



l/3r- 



-2/3 



lO'-'cm" 



^6 'T^(Azt,)-^'^ (52) 

'^^ for our typical cloud in order 
to match observations, this would imply nici ~ 3 x 10^ Mq. 
The angular size of the cloud is related to the cloud 
radius by 



10" 



D 



deg 



(53) 



where D is the distance to the cloud. If we set D — Re 
feauation ll8|l . then the typical cloud area on the sky will be 

6,1-0.25 deg' m^/'^-'''^ (54) 

In order to match a typical size of dd — 0.5deg' we will 
need rrici ~ 3 x 10^ M0, which is nicely in line with what 
we needed to match the column density above. Of course, 
the angular sizes observed correspond to HI sizes, and one 
might expect the outer radius in neutral material to be 
some what smaller that the full cloud radius. In the mod- 
els of lMalonev fc PutmarJ ll2003ll . rni — 0.7rci for constant- 
density clouds similar to the type we consider here. If we 
adopt this assumption, then the coefficient in equation 1541 
would scale to 0.18 deg', and push our preferred mass to 
mci ~ 5 X 10*^ Mq. 

Finally, there are roughly 2000 HVCs in the HIPASS 
sample covering the southern sky. If we double this, we 
can estimate that the full halo should contain ~ 4000 such 
clouds. The number of clouds we expect in the halo is simply 
A'ci — Mci/rrici. In the previous section, we assumed cloud 
masses of md = 5 x 10® M© and computed that the total 
residual cloud mass in the halo would be Aid — 2 x 10^" M©. 
This implies Nd = 4000, consistent with the number ex- 
pected from the HIPASS count. 

It is remarkable that the numbers in all of these cases 
work out to favor roughly the same cloud mass, rud — 
(3 — 5) x 10^ Mq. This is likely something of a coincidence 
considering the crudeness of our model. Although we have 
not focused on it here, the HVCs are observed to have a 
distribution of sizes and column densities. This might be 
achieved by allowing a distribution of cloud masses and 
some more sophisticated treatment of how they might dis- 
rupt upon approaching the galaxy. Nonetheless, we take it 
as a positive sign that our simple model is able to match the 
rough characteristics of HVCs using a single cloud mass. 



8 QUASAR ABSORPTION SYSTEMS 

It has long been assumed that quasar absorption systems 
can be identified with the gaseous content of galaxy halos 
JBahcall fc Spitzer 1969) . High column density systems like 
Lyman limit and CIV systems ar e observed to have nearby 
optical counterpart s f or z < 1 l|B e rg eron &: Boisse| |l99l|: 
[ Chen et al.l l2001allbl: ISteidel et alj Il997l: iLanzetta et alJ 
Il995r . Theoretically, these systems have been modeled as 
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Figure 10. The covering factor to CIV absorbers as a function 
of impact parameter for our fiducial model galaxy at 2 = 0.5 
(solid line). The er ror bars are representative of the results of 
IChen et alj 120013^ . The sharp truncation that is indicated by 
the data is roughly at R^ which is what would be expected in our 
model. More detailed comparisons are difficult because the data 
spans a wide range of redshifts and luminosities. However, with 
more data absorption systems will provide important constraints 
on the properties of the warm clouds in different mass halos. 



arising from warm clouds emb edded within a hot galaxy halo 
||Mo fc Miralda-Esc udelll996ll in a way that is quite similar 
to what wc desc ribe he re. 

IChen et alJ (120013^ find in their sample that when the 
impact parameter of the quasar is < 97 kpc, 67% of galax- 
ies show CIV absorption, while when the impact parame- 
ter is > 97 kpc, only 6% of galaxies show CIV absorption 
systems. We compare this rough expectation to the cloud 
covering factor as a function of radius calculated using our 
fiducial parameters for a galaxy at z = 0.5 (Figure tTUIl . In 
our model there is also a sharp drop in covering factor that 
occurs at the cooling radius of the halo. Thus our model is 
in qualitative agreement with the observations. 

To properly model quasar absorption systems one must 
be able to connect galaxy luminosity and type to a halo's 
maximum circular velocity as a function of redshift. Then 
each observed galaxy's gaseous halo can be modeled and 
compared to observations. With a great deal more absorp- 
tion data it will be possible to constrain the masses and 
numbers of clouds in a halo as a function of halo mass and 
redshift. One quantity that would be useful to know is the 
average relationship between the column density of the ab- 
sorption system and the total amount of mass along the line 
of sight. This should be possible combining weak gravita- 
tional lensing with a largo sample of quasars and absorption 
systems (Mailer et al. 2002), and some progress has been 
made on this front iJMenard fc Perouxl l2003i l . We discuss 
how other properties of absorption systems may be useful in 
constraining the properties of the warm clouds in HIOI 



9 THE LUMINOSITY FUNCTION 

One of the fundamental goals in galaxy formation model- 
ing is to understand why there are so few galaxies with 
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baryonic masses larRer than ~ lO^^^M© . The coohng-time 
arguments (e.g. IWhite fc ReesI Il978l) were originaUy de- 
vised with the goal of explaining this upper-mass cutoff, 
but it is now accepted that the coohng radius treatment 
alone cannot do the job in the contex t of modern ACDM 
cosraoloKv (e.g. LffiMte , & Frenk] ll^gl): [Thc^r^l fc ' ^einbere 
I1995I: ISomerville fc Prim ack 19991: iBenson et alJl2003l) . The 
ACDM halo mass function (or velocity function) follows a 
near power-law distribution over t he mass-scale (or velo city 
scale) of the Milky- Way halo (e.g. lGonzalez et alJlJOOcl) . In 
contrast, the luminosity function drops quickly above the lu- 
minosity scale of the Milky Way (e.g. Blanton ct al. 2003) ■ 
The cooling radius treatment reduces the fraction of gas that 
cools in high-mass halos, but only moderately, and certainly 
not at the level required to explain the characteristic lumi- 
nosity of galaxies (see Fig. 1121 below). 

The issue i s highlighted noticeably by the results of 
iBell et al.l (|20o3) who used data from 2MASS and the SDSS 
to construct the baryonic mass function of stars-|-gas in 
the local universe. They concluded that the number den- 
sity of galaxies falls off sharply above a cold baryonic mass 
of M* ~ 1O"M0 (shaded band in Fig. EU. By integrating 
their mass function, they found that the total mass in cold 
baryons in the local universe is only ~ 10% of the total bary- 
onic mass expected from BBN and the concordance ACDM 
model. While the original cooling arguments suggested that 
most of the baryonic mass would end up in stars, it seems 
now that most of the baryons have ended up in hot gas, or at 
least in some state that is no t associated with central galax- 
ies. As discussed bv lBenson et al. (2003), explaining the the 
sharp cutoff at the bright end of the luminosity function is 
difficult within the standard scenario without resorting to 
extreme conduction (above the Spitzer value) or hot super- 
winds with energies beyond expectation. 

Figure HTl shows how this problem might be alleviated 
by allowing a two-phase medium to develop during the cool- 
ing process. The dashed line shows the total baryonic mass 
within the cooling radius, M_r^ , as a function of halo maxi- 
mum circular velocity. This mass, associated with the central 
galaxy in the standard treatment, continues to rise rapidly as 
the halo velocity increases and would naively lead to a popu- 
lation of giant ~ 8x 10^^ Mq galaxies associated with galaxy- 
group halos with Vmax — 500km s~^. The dot-dashed line, 
showing the total cooled mass, Mc = Mr^ — Mh, determined 
by the two-phase treatment described in 2] is more encour- 
aging. Interestingly, this mass approaches a characteristic 
value of Mc ~ 2 X lO" M© for halos with Knax^ 250km s"^ 

The reason why the total cooled mass approaches a con- 
stant in large halos is that the hot halo core, Afh, grows 
rapidly with Vmax because of the corresponding increase in 
the cooling density (short dashed line). This compensates for 
the increase in Mr^ , resulting in Mc = Mr^ — i\fh -^ con- 
stant. Once Knax^ 250km s~^, -Re ^constant (Fig. I^J and 
Mh ~ pcRc oc pc oc Knax (see equations 1181 and 1121 and Fig. 
|5J|. In this regime, the total mass inside the cooling radius 
for our initial hot gas profile also increases proportionally 
to Knax, so that Mr^ OC Vmax as wcU. Therefore, both the 
total mass within Re and the mass in the hot core increase 
with Vmax in the same way (the dashed lines in Fig. IllH . 
The amount of cooled mass remains constant at the value it 
had when the slopes began to match. Of course for different 
assumptions about the initial hot gas profile the slopes may 
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Figure 11. The total baryonic mass within the cooling radius, 
Mr^ (long-dashed line) as a function of halo Knax , assuming a gas 
metallicity oi Zg =0.1. Also shown is total mass that we expect 
to have cooled, Mc (dot-dashed line) and the total mass in the 
hot core (short dashed line). The total cooled mass approaches a 
constant in high- Vmax halos because much of the mass within the 
cooling radius ends up in a pressure-supported hot-gas core. The 
mass in cooled gas that accumulates into the central galaxy Mg is 
plotted for two different assumptions about the cloud masses. The 
solid line shows the resultant galaxy mass calculated assuming a 
constant md = 5 X IO^^Mq and the dotted line shows the result 
assuming that the cloud mass is set by conduction (see text). 



not exactly match and the amount of cooled mass may in- 
crease (or decrease) slightly instead of remaining constant. 
However, it still will change much less drastically than the 
total mass within the cooling radius. 

The solid line in Figure 1111 shows our expectation for 
the central galaxy mass as a function of Vnax. This is de- 
termined using the methods outlined in H6.2I assuming a 
constant cloud mass of nici = 5 x 10^ Mq. We see that 
in this case, the total mass in cold gas that ends up in 
the central galaxy approaches a value Mg ~ 1.5 x 10^^ M© 
for halos with Vmax^ 250km s~^. Of course, this treatment 
has not included any mergers between halos, so that this 
maximum galaxy mass really represents the max;imum mass 
galaxy sitting within a relatively "quiescent" halo. An in- 
teresting implication is that forming a galaxy more massive 
than ~ 1.5 x lO^^M© would require a merger. This may 
be relevant in explaining why spheroidal galaxies tend to 
dominate the bright-end of the luminosity function. 

Another possibility is that the characteristic cloud mass 
is not constant, but scales with the halo temperature in some 
way. The short-dashed line assumes that cloud masses are 
set by conduction (md oc T^^''"') as described in >I5.2I (" dot- 
dashed line in Fig. |SJ|. In this case, cloud masses become 
quite large, rrid^ 10* M© in Vmax^ 300km s~^ halos. The 
infall time of clouds scales as rin oc mj (equation 1431 1441 
and thus becomes quite long at T increases. Clouds tend to 
remain in the hot gas halo rather than fall in to contribute to 
the central galaxy in this case. Interestingly, massive clouds 
of this type will likely be Jeans unstable in high-mass halos 
(e.g. Fig. ISJ. In this case, they may form stars and become 
small galaxies on their own (see our discussion in HIOI I. As 
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Figure 12. The cu mulative baryonic mass function of galaxies 
reported by iBell et al. 1 2003) (shaded band) compared to the 
cumulative mass function of halo baryons (top solid line). The 
short-dashed line shows the (central) galaxy mass function that 
arises from assuming all of the mass within each halo's cooling 
radius cools onto the central galaxy, Mn^ . The dot-dashed line is 
the cooled mass (Mc) function that arises in our picture, which 
allows for the presence of a hot corona in each halo (see text and 
Fig, nil . Finally, the lowest solid line shows the central galaxy 
mass function that results from modeling the survival probability 
of cooled clouds in the halo, assuming a typical cloud mass of 
5 X 10^ Mq. Only clouds that fall to the center of each halo are 
assumed to contribute to the central galaxy. As mentioned in the 
text, no merging has been accounted for in this estimate. Merging 
will tend to populate the massive tail of the galaxy mass function, 
likely bringing it even more closely in line with what is observed. 



in the fixed rrici case, giant CD galaxies could only form via 
mergers in this picture. 

The mass function of galaxies is shown in Figure^] The 
shaded band shows the baryo nic mass function of galaxies 
in the Universe determined bv lBell et alJ i2003r . The width 
of the band indicates their uncertainty (which comes mainly 
from the IMF). Compare this to t he upper solid line , which 
shows the halo mass function of ISheth fc TormenI (119991) 
scaled by the mass of baryons in each halo (Mv — > fbM^). 
The offset is roughly a factor of -^ 10 in normalization, and 
from this one can immediately see most of the baryons in 
the universe cannot be in the form of cold, galactic mate- 
rial. The mass function that results from assuming that all 
of the baryons within each halo's cooling radius cool onto 
a galaxy (shown by the short-dashed line) cannot solve the 
problem. That is, it shows no sharp drop in galaxy counts 
above ~ 10^^ Mq. As expected from the above discussion, 
the "cooled mass" function derived using our multi-phase 
picture does much better in accounting for this cutoff (dot- 
dashed line). 

The mass function of gas that we expect to actually fall 
into the central galaxy is shown by the lowest solid line in 
Figure^] In this estimate we have assumed that clouds have 
a typical mass of 5 x 10^ Mq (as in the solid line in Fig. lllll . 
The high-mass cutoff compares quite well to the data in this 
case. We expect galaxies more massive than this cutoff to be 



produced solely via mergers, and we suggest that mergers 
will tend to populate the tail of the mass function above 
~ 10^^ Mq, bringing it even more in line with observations. 
As is clearly seen, our cooling scenario will not help explain 
the well-known faint-end slope problem (the low number 
density of galaxies smaller than ~ 5 x lO^^M©). Of course, 
feedback likely plays a major role in this regime. 

It is quite clear from Fig. 1111 that the galaxy mass 
stops increasing because most of the baryons remain in the 
hot core. The hot gas in higher temperature, more massive 
halos can be observed with x-ray telescopes. These obser- 
vations suggest a more complicated picture then we have 
been describing her e; where energy injection from AGN (e.g. 
lOmma et alJl2004l) or other sources may be needed to ex- 
plain the obse rved x-ray luminosity vrs. temperature rela - 
tionship (e.g. iMushotzkv fc Scharl ITqQtI: IWu et alJlJOOol) . 
Energy injection into the hot gas may also occur in galaxy 
mass halos. We suggest that it may not be needed to explain 
the high mass cutoff in the luminosity function, although 
undoubtedly heating processes occur to some extent. 



10 IMPLICATIONS, OBSERVATIONS, AND 
FUTURE WORK 

Central to our model is the existence of a hot, low den- 
sity med ium tha t surro unds galaxies — an idea first pro- 
posed bv lSpitzeil 1119621) . The hot gas density we expect for 
a Milky- Way type system is Uh ~ Ph/^irrip ~ 8 x 10~^cm~'^ 
at ~ lOOkpc. Evidence that such a corona exists has been 
growing in recent years. For example, gas clouds in the Mag- 
ellanic Stream are more easily understood if they are c on- 
fined by a hot gas medium JStanimirovic et al.l |2002^ as 
are the shapes of supergiant shells along the outer edge 
of the LMC ^e Boer et al. 1998). Indeed, detection of H- 
alpha emission at the leading edges of clouds in the Mag- 
ellanic Stream is best explained by ram-pressure heating 
from a substantial hot gas halo (n^ ~ 10~*cm~"^) at ~ 
50kpc ('Weine r fc Williama Il99(i) . Unfortunately, most of 
the quantitative limits on galactic hot gas densities are 
model dependent, but the currently limits are at least in 
line with our predictions: Uh^ (1 — 10) x 10~^ cm~'^. (e.g. 
Snowden eral]ll997l: [Benson et ai]lj200(]l:JBiitz fc Robishawl 



ISnow 
|2000: 



Moore & Davis 1994; Murali 2000; "Tr ipp et al.ll2003l) 
Interestingly, Quilis fc Moore (2001) argue that gas densi- 
ties as low as ~ 10~*cm~^ cannot pr oduce the head-tail 
posit ion- velocity gradients observed by iBriins et al.l (l2000l . 
I2OOII) for ~ 20% of HVCs. Our hot halo core is somewhat 
denser than this, so the head-tail gradients may be consis- 
tent our scenario. 

One of the most promising methods for probing the hot 
gas halo is to observe it in absorption. Indeed, OVI in or 
around th e Milky Way halo has been detected in this manner 
by FUSE dSavaee et alJbOOd; ISembach et al.ll2003D . Higher 
ionization lines (e.g. OVII a nd OVIII), can be ob served by 
XMM-Newton and Chandra. iNicastro etlii] (l2002h reported 
the first detection of highly ionized oxygen and neon in prox- 
imity to the Local Group and similar detections have fol- 
lowed. 

There is some debate over exactly where this absorp- 
tion takes place and how the obs erved OVI and OVII ab- 
sorbers are related (e.g. compare ISembachI 120031 : iNicastrd 
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l2003f) . In the context of the corona we predict, measure- 
ments to known sources in the Local Group may help avoid 
confusion. For example, our fiducial model predicts a col- 
umn density in Hydrogen of Nh = lO^^cm"^ along a line 
of sight to the LMC, assuming a distance of 50kpc. With 
Zg = 0.1, this gives a total expected column density in oxy- 
gen of 5.1 X lO^'^cm"^. At the temperature and density we 
expect for the Milky Way corona, most of the oxygen should 
be in the form of OVII (e.g. Mathur et al. 2003' ). but observ- 
ing multiple elements and ionization lines would be useful for 
constraining the precise ionization state, temperature, and 
density of the hot medium(s) responsible for any absorption. 

A related test will come from searches f or metal lines 
associated with HVCs. ISembach et al.l i2003l) have argued 
that high-velocity OVI features observed by FUSE high- 
light the boundaries between warm clouds of gas and a 
highly extended, hot, low-density corona around the Galaxy. 
A related analysi s suggests that th ese systems are associ- 
ated with HVCs JTripD et al.ll2003l) . We would expect just 
this situation in our model. Note that lNicastro et alJ (120031) 
have argued that the high-velocity OVI absorbers are better- 
described by a Local Group population, but they also al- 
low for the possibility that they trace an extended Galactic 
corona. This second interpretation is in line with our expec- 
tations. 

The soft x-ray background provides a less-direct method 
for probing the hot gas cores of galaxy halos. In our model, 
the hot-gas cores are expected to be of low density, and 
to have a rather low x-ray surface brightness, not directly 
detectable by current x-ray satellites. However, this hot gas 
will contribute to the soft-x-ray background. Predictions for 
the contribution to the soft-x-ray background may provide 
interesting limits on the model. 

More detailed comparisons with existing and future 
HVC data will require more realistic models of clouds in 
a hot halo. The properties of clouds will need to be modeled 
in the presence of an ionizing background field and interac- 
tions with the hot gas bac kground should be p roperly taken 
into account. For example. IWeiner et alJ |(200?) have argued 
that measurements of H-alpha emission in HVCs can put 
constraints on their distances as long as H-alpha recombi- 
nation is caused by photoionizing radiation from the Milky 
Way. Another possibility is that the H-alpha recombination 
is due to coUisional ionization caused by ram pressure in- 
teractions with the hot gas halo. Such a scenario would be 
consistent with the OVI observations discussed above. Fur- 
ther, the ambient pressure from the hot gas halo is expected 
to vary as a function of radius from the halo center, and 
this would lead to varying cloud sizes (and densities) at 
fixed mass. The spectrum of HVC sizes and column den- 
sities might then be used to constrain the nature of the hot 
gas corona, and even to test the mass spectrum of clouds. Of 
course, without knowing the distances to individual HVCs, 
this can only be done in a statistical sense. 

Searches for clouds around other galaxies will help es- 
tablish whether clouds of the type we discuss are as com- 
mon and numerous as we expect, and might even be used 
to test how cloud masses and hot c ore properties vary with 
halo mass and galaxy luminosity. IPisano et al.l i2004) re- 
cently performed a search for HI clouds around three nearby 
galaxy groups and found that if a population of HI clouds 
exists around galaxies like the Milky Way, they must be 



clustered within 160 kpc and have HI masses ^ 4 x W^ Mq. 
The clouds we expect are consistent with these limits, but 
should be detectable if the detection limits are relaxed only 
slightly. Thilker et al. ( 2004) have discovered a population 
of ~ 50 HI clouds around M31, with HI masses of ~ 10*^ Mq. 
Our model would suggest that these are likely the most mas- 
sive of the many thousands of clouds that should surround 
M31. We predict that deeper surveys with better angular 
resolution will find these clouds. 

Quasar absorption systems provide another important 
avenue for determining the properties of warm clouds. While 
cloud populations of this type cannot be studied in detail, 
the study of absorption systems can probe wide range of halo 
types. A cross-correlation between absorbers and galaxies 
may yield useful information on cloud sizes, densities, and 
covering factors. In the future, a cross-correlation survey, 
similar to that of Chen et al. (2001a), but utilizing a large 
optical survey, would yield tight constraints on the distribu- 
tion and scaling of cloud masses. 

Another advantage of absorption systems is that they 
probe the gaseous halos of galaxies at early stages of 
formation, z ^ 3. Multi-phase cooling may be a cru- 
cial ingr edient in und e rstand ing the properties of these 
systems. iMaller et al.l (|2003r pointed out that a large 
fraction of the halo gas must be in form of warm 
clouds in order to explain the observed kinematics of 
the h igh-ion component in d amped Lyman alpha sys- 
tems dWolfe fc Prochaskall200(ll . Including the multi-phase 
medium self-consistently will be important for precise com- 
parisons with this data. Indeed , damped systems them- 
selves show complex kinematics JProchaska fc Wolfelll997l . 
Il998f) that cannot be explained in CDM cosmologies with- 
out the presence of a l arge amount of gaseous substructure 
HHaehne lt et alj Il99^ : iMcDonald fc Miralda-Escud J IT999I: 
Mailer et al. 2001J). Warm clouds may make an important 
contribution to this substructure, possibly after they are dis- 
rupted by tidal forces. 

If pressure-supported clouds exist in the halos of most 
galaxies, this could be important for interpreting the flux 
ratios of multiply-image d quasars (jMetcalf fc Madau 200J| ; 
iDalal fc Kochanekl 12003 : iMoustakas fc Metcal J l2003ll . The 
flux ratio anomalies have been used to argue the existence of 
the low-mass dark matter halos predicted by ACDM N-body 
simulations dMoore et al.lll999L IkIvpIu et alJll999ll . How- 
ever, there are some indications that the fraction of mass 
in low-mas s ^ lO^Mp substructures is even higher than 
predi cted dZentner fc Bullockl 120031: IMoustakas fc Metcall 
l2003h . The warm clouds we have described here will also 
cause fluctuations in the gravitational potential, and may 
be important. Indeed, the mass fraction in clouds is ex- 
pected to be at the few percent level, and this is roughly 
the same as the dark matter substructure population. Fortu- 
nately, the warm clouds may also be detected by absorption 
in the quasar spectra so it may be possible to disentangle 
the two signals. Also, since the clouds are expected to be 
roughly constant density, they may not provide as strong a 
signal as the more concentrated dark matter clumps. 

Multi-phase cooling might also help resolve the long- 
standing problem of forming disk galaxies without an- 
gular momentum loss in cosmological simulations (e.g 
iNavarro fc Steinmetd 120001) . Specifically, if cooled gas re- 
mains in warm clouds instead of settling into the galaxy. 
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then those clouds can retain and gain angular momen- 
tum during mergers. When t he clouds eventually fall in, 
they will produce large disks (Mailer & Deker2002") . This 
scenario would be especially helpful if angular momen- 
tum in dark matter halos is predominately a cquired in 
merge rs llMaller et alJl2002l: IVitvitska et al.ll20n2h . Interest- 
ingly, iRobertson et al.1 i2004f) showed that by allowing a 
cold/warm medium to exist within cooled, star-forming ma- 
terial, they could improve the likelihood of disk formation in 
cosmological simulations. Specifically, the disk is more stable 
to its own self gravity, and less likely to fragment and loose 
angular momentum after it forms. Of course, this effect will 
only help if the material that forms the disk initially retains 
a large amount of angular momentum. It is in the retention 
of halo angular momentum that the warm/hot cloud picture 
becomes important. Therefore, a full multi-phase approach, 
with an allowance of both cold/warm and warm/hot phases, 
could lead to more success in this direction. However, as 
we mention in Hill there are significant computational chal- 
lenges to overcome. 

Finally, as seen in Fig. |S1 the Jeans mass for a cloud 
decreases as a function of halo temperature, making it more 
likely that a cloud will collapse under its own gravity in high- 
temperature halos. This possibility will be more likely if the 
typical cloud mass increases in high-temperature halos. This 
is what is expected, for example, if conduction sets the char- 
acteristic cloud mass as discussed in US. 21 If a cloud's mass 
exceeds the Jeans mass then it will likely fragment to form 
stars. It is perhaps to be expected then that there should 
exist a population of low-mass galaxies, born of fragmented 
gas in clusters, with no associated dark matter. Dwarf galax- 
ies of this type would likely be younger than other low-mass 
spheroids, and have relatively low mass to light ratios by 
comparison. 



11 CONCLUSIONS 

In this paper we have taken a step towards modeling the 
complex realities of astrophysical hydrodynamics using a 
simple analytic treatment that allows the development of a 
two-phase warm/hot medium during gas cooling. Appealing 
to standard cooling instability arguments, we showed that if 
cooling proceeds by the formation of warm clouds embedded 
within a low-density hot gas background then this can ex- 
plain the characteristic upper limit in the observed baryonic 
masses of galaxies, ~ W^^ Mq (©. In the standard treat- 
ment, all of the mass within the cooling radius of each halo 
cools onto the central galaxy, while our approach allows the 
survival of a hot gas core with a density close to the cooling 
density in each the halo. The fraction of mass that remains 
in the hot core component is large in high-mass halos be- 
cause the cooling density is high, and this gives rise to an 
upper-mass limit in cooled material in these systems. 

When applied to Milky- Way size halos, the standard 
single-phase treatment over-predicts the Milky Way mass 
by more than a factor of two, while our multi-phase treat- 
ment helps explain the Milky Way galaxy mass naturally, 
without the need for excessive feedback ([EJ. Because of the 
thermal instability, we argue that galaxy formation should 
proceed via the infall of warm, pressure-supported clouds. 
Further, if the typical cloud mass is nici ~ 5 x W^Mq, the 



residual cloud population is significant, and we identify these 
pressure-supported fragments with the observed High Ve- 
locity Cloud population of the Milky Way (0. The typical 
Galacto-centric distance to HVCs in our picture is set by 
the cooling radius, ~ lOOkpc. The same cloud mass helps 
explain the baryonic mass of the Milky Way, can account 
for high- ion absorption systems in distant galaxies {^^, and 
allows clouds to survive destructive processes in the halo 

Including the multi-phase treatment presented here in 
standard semi-analytic models should be straightforward. 
Currently this type of modeling tracks two phases of gas: 
the mass contained in hot halo gas, and the mass in the 
"cold" central galaxy. To include the multi-phase medium, 
one must simply add the hot core from eauation l24l to the gas 
in the hot phase, and include an additional accounting for 
warm cloud material. As sketched in equation 1461 gas first 
cools into warm clouds and warm clouds become deposited 
in the central galaxy on an infall timescale. Once clouds fall 
in, the typical recipes for star formation may be applied, 
although it is likely that less feedback will be needed. 

Comparison with cosmological hydrodynamic simula- 
tions may be somewhat more difficult. At present, these sim- 
ulations do not resolve the multi-phase structure of halo gas, 
and this may lead them to predict infall rates of cooled gas 
that are similar to those expected from simple cooling radius 
arguments. In order to test these expectations, simulations 
would need to resolve typical cloud masses of ~ 10^ Mq as 
well as properly follow clou d fragme ntation in the diffuse 
gas ha lo. I Yep es et al.l ilQQTT and later ISpringel fc Hernouisli 
(2003) have in fact included a sub-grid model for multi- 
phase gas in their SPH codes, but this applies only to the 
star-forming, cold/warm medium. This approach seems to 
alleviate many of the problems faced by similar codes in 
the past (e.g. IRobertson et al.l 120 04'). but not the "over- 
cooling" problem discussed here. If the computational chal- 
lenges can be overcome, hydrodynamic simulations with full 
multi-phase cooling may yield even more encouraging re- 
sults. 

In conclusion, multi-phase cooling is expected on theo- 
retical grounds and can alleviate many of the problems that 
arise in the standard, single-phase procedure. The fact that 
warm clouds of the type predicted seem to be observed only 
further enforces the relevance of adopting this approach in 
models of galaxy formation. If we are correct, then the High- 
Velocity Clouds of the Milky Way are tracers of the funda- 
mental fuel supply that governs galaxy formation in the Uni- 
verse. Set in this context, the study of hot gas, HVCs, and 
their counterparts in external galaxies will have significant 
impact on how we understand galaxies and their assembly. 
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Table Al. Properties of the cooling function for several different 
choices of hot gas metallicity, Zg , as parameterized by the fitting 
function in equation IA2I The temperature T;, is the temperature 
above which Bremsstrahlung radiation begins to dominate cool- 
ing, and Vmax(Ti,) is the corresponding halo maximum circular 
velocity. The parameter a represents the fitted slope to the cool- 
ing function in the regime: Tr <T <Tm, where Tr = 1.5 X 10''K 
and Tm ^ 1.5 X lO^K for Zg > and formally T^ = T^ for 
Zg = 0. Finally, Az is the value of the cooling function in units 
of 2.6 X 10~^^cm'^ergs^^ as defined in equation 1131 
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APPENDIX A: A. THE COOLING FUNCTION 

Astrophysical plasmas with temperatures greater than W^K 
primarily cool by radiative processes. The cooling funct ion 
A(r) can be calculated (e.g. ISutherland fc Dopitalll993i as 
a function of the gas metallicity Zg . It is useful to introduce 
the dimensionless cooling function L23(T): 
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A(r) = 10"^^cm='ergs-^L23(r) (Al) 

Figure0shows A(T, Zg) plotted as a function of temperature 
for five example metallicities. Note that for galaxy-sized ha- 
los (Vmax — 100 — 200km s~^) and mildly-enriched gas, the 
dimensionless cooling function takes values of order unity. 

We can approximate the cooling function by a series of 
power-law fitting functions that captures the important scal- 
ings. For zero metallicity gas there are two important tem- 
perature regimes: T > Th, when the dominant cooling pro- 
cess is Bremsstrahlung radiation, and Ti, > T > Tr — 1.5 x 
10*K, where the dominant cooling process is the recombina- 
tion of hydrogen. For enriched gas, there is a third important 
temperature scale, Tm, where metal line cooling becomes im- 
portant. The Bremsstrahlung region becomes important at 
higher temperatures for more metal rich gas, and we find 
that to good approximation Tt = 10^ + 1.5 x lO^Zg^/^K. We 
also adopt T™ = 1.5 x lO^K. 

The cooling curve is then given by, 

L23-^(r) = -- ^T 



12 






Tr<T<T„ 



L23^ (T) = L23-' (T, 



L23-^(r) = L23-f(r5)(|r 



T,n <T<Tt 



1/3 



T>Tt, 



(A2) 



where a 



ilnZ, 



[Zg 7^ 0), and we have added the 
superscript / to indicate that this is a fit to the true cooling 
function. In the case of zero metallicity gas a — —0.8 and the 
middle expression is not used (that is Tm ~ Tt). Values of a 
and T(, for some example metallicities can be found in Table 
lAll This approximation is at worst good to within a factor 
of 3 and is much better than that for most metallicities and 
temperatures. 



APPENDIX B: B. HYDROSTATIC 
EQUILIBRIUM SOLUTIONS 

Here we derive the hydrostatic equilibrium solutions for gas 
profiles, assuming that the gravitational potential of the sys- 
tem is dominated by an NFW background halo. The hydro- 
static force balance equation is 

dP _ ~V^(R)p{R) 
dR 



R 



(Bl) 



where R is the distance from the center of the spherically- 
symmetric halo. The density follows that given by equation 
|S] and the implied rotation curve follows 



V^{R) = 



GM{R) 



9.26c; 



Rsf{R/Rs 



(B2) 



R " R 

where f{x) = ln(l + x) ~ x/{l + x), and we adopt Cg = 
Vmax/\/2 as the sound speed written in terms of the maxi- 
mum circular velocity. 

If we assume P — K-^p^' then K^ = Cgp].~~' and the 
solution to the hydrostatic equation for 7 7^ 1 is 

9.26(7-1) /ln(l + x-) ln(l + Cc 



Pg{x) 



Pc 



1 + 



7 



Cc 



For 7 = 1 the solution is 

pg{x) = PcCXp 



9.26 , ,, , 9.26 , ,, ^ 

ln(l + x)- -—- In 1 + Cc 

X Gc 



(B4) 



where Cc = Rc/Rs- In writing this solution, we have re- 
quired that the hot gas density at the cooling radius, Re, is 
equal to the cooling density, pc- For the adiabatic assump- 
tion adopted in the paper, 7 = 5/3 and the solution for 
pressure, density, and temperature becomes 



PA^) -- 


= Pc 


Pa{x) = 


= Pc 


T,{.x) -- 


= r 



3.7 



ln(l + 2;) 



1 + — ln(l + x) 

X 



3^ 

Cc 

3^ 

Cc 



ln(l + Cc) 



ln(l + Cc) 



5/2 



,(B5) 



3/2 



1 + — hi(l + x) - ^ ln(l + Cc 

X Cc 

For our adopted adiabatic profile, the pressure, temper- 
ature, and density increase slowly towards the center of the 
halo, reaching core values as x — > 0: To = aT, po — a'^'^pc, 
and Po = a^^^Pc with q = 4.7 - 3.7Cc"^ ln(l + Cc). For the 
typical range Cc = 2 — 20 we find a rather modest range of 
values a ~ 2.7 — 4.1. It is worth pointing out that while the 
central density is higher than what we have defined as the 
"cooling density" (set at T = Tc), it does not imply that the 
system is drastically unstable to cooling. Indeed, the cen- 
tral hot gas is likely to have reached its state of density and 
pressure more recently than the non-radiating gas at large 
radii. Moreover, for Tm < T < Th, the local cooling time 
will scale as TcooI oc T/pA{T) oc T^/p oc a^'^ , and thus the 
central gas will cool more slowly than the outer halo gas. 

Most of the gas, by volume, is quite close to the state 
at Re- The total mass in the form of hot gas takes the form 



Mb = Yri?pcr?d(Cc) 



(B6) 



where the function i]d is determined by numerical integration 
to be well fitted by 



r,d(c) ~ 1.42c°-^ [1 + (c/1.65)i-'] 



(B7) 



which is good to < 1% for c — 1 — 20, and is a rather weak 
function of Cc: r/d ^ 1-45 - 1.25 for Cc = 1 - 20. Clearly the 
average gas density within re will be ph = rjdPc- In the main 
part of the paper we work under the approximation that rjd 
is a constant and adopt a typical value of ryd = 1.35. 

Similarly, we can estimate the volume-averaged pressure 
of the gas within Re : 

47r f^" PMR^dR 
a = -\^^^37^=P...(Cc), (B8) 

where the function rjp is found to be well-fitted by 
r;p~2.5c°-^^[l + (c/3.9)^-«]-°-'", (B9) 

and spans the range rjp ~ 2.1 — 3.3 (for Cc = 1 — 20). As with 
the gas density, in the main part of the paper we work under 
the approximation that the hot gas pressure can be well- 
represented by a constant, with Ph = rjpPc, and adopt the 
typical value r]p ~ 2.7. The temperature profile is relatively 
fiat, and we adopt Th = T throughout the paper. 
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